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Abstract 

Atomistic  simulation  with  a  modified  version  of  the 
DYNAMO  code  was  used  to  calculate  the  point  defect  proper¬ 
ties  of  the  B2-type  intermetallic  compounds  NiTi  and  FeTi. 
The  calculated  energies  are  believed  to  be  experimentally 
accurate  to  within  +/-  0.1  eV. 

The  antisite  pair  formation  energies  for  the  alloys 
were  calculated  to  be  1.015  eV  for  NiTi  and  0.726  eV  for 
FeTi.  The  antisite  defect  and  vacant  Ni  or  Fe  lattice 
sites  were  found  to  be  weakly  bound  with  a  binding  energy 
0.094  eV  for  NiTi  and  0.055  eV  for  FeTi. 

The  vacancy  defect  configurations  in  NiTi  and  FeTi  were 
idcr.t i  :_1 .  Vacant  Ni  or  Fe  lattice  sites  were  preferred 
and  had  energy  values  of  6.356  eV  for  NiTi  and  5.899  eV  for 
FeTi.  The  Ni  and  Fe  vacancy  migrations  had  energy  barriers 
of  1.272  eV  for  NiTi  and  1.737  eV  for  FeTi.  Fa.-ovir.g  a  Ti 
atom  resulted  in  a  neighboring  antisite  defect  caused  by 
migrating  Ni  or  Fe  atoms  into  the  vacant  Ti  lattice  site. 
These  configurations  had  energy  values  of  7.024  eV  for  NiTi 
and  6.336  eV  for  FeTi.  The  vacancy  defect  formation  ener¬ 
gies  were  calculated  to  be  1.48  eV  for  NiTi  and  1.07  eV  for 
FeTi . 
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The  interstitial  defect  configurations  consisted  of 
Ni-Ni  or  Fe-Fe  split-interstitial  dumbbells  centered  on  a 
Ti  site  with  one  or  two  adjacent  Ti  antisite  defects.  Ti 
interstitials  resulted  in  two  adjacent  Ti  antisite  defects 
with  a  < 0 .1 1  >  direction  Fe-Fe  dumbbell  or  a  <111>  direction 
Ni-Ni  dumbbell.  The  Ti  interstitial  defect  energies  were 
-2.395  eV  for  NiTi  and  -1.558  eV  for  FeTi.  Ni  or  Fe 
interstitials  both  resulted  in  a  single  adjacent  Ti  anti¬ 
site  defect  with  Ni-Ni  or  Fe-Fe  dumbbells  oriented  in  the 
<111>  direction.  The  Ni  or  Fe  interstitial  defect  energies 
were  -2.602  eV  for  NiTi  and  -1.945  eV  for  FeTi. 
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CALCULATED  POINT  DEFECT  PROPERTIES  OF  THE  INTERMETALLIC 
COMPOUNDS  NICKEL  TITANIUM  (NiTi)  AND  IRON  TITANIUM  (FeTi) 


I .  Introduction 


Background 

The  study  of  point  defects  for  several  decades  has  been 
of  interest  due  to  the  role  defects  play  in  determining  the 
physical  properties  of  crystalline  solids.  This  is  espe¬ 
cially  true  in  those  defects  controlling  the  transport  of 
matter  and  the  properties  stemming  from  it  (Crawford,  1972, 
1) .  For  example,  atom  vacancies  in  the  lattice  structure 
are  accepted  as  the  dominant  mechanism  for  self-diffusion 
in  metals  (Agullo-Lupez ,  1988,  196).  A  point  defect  can  be 
defined  as  "a  perturbation  that  moves  as  an  entity  and 
whose  properties  are  always  the  same,  independent  of  how  it 
has  been  generated  and  where  it  is  located  in  an  otherwise 
perfect  crystal"  (Gruber,  1964,  38).  Some  of  the  first 
papers  in  the  study  of  point  defect  properties  were  done  in 
calculating  the  defect  formation  and  migration  energies  by 
Huntington  and  Seitz  in  1942  (Gruber,  1966,  21).  Irradi¬ 
ation  of  solids  with  electrons,  neutrons  and  heavy  ions 
result  in  lattice  point  defects.  Thus,  experiments  with 
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irradiated  materials  have  been  used  in  the  study  of  point 
defect  formation  and  the  development  of  atomistic  simu¬ 
lation  models  of  point  defects  (Gruber,  1966,  21). 

Recently,  theories  in  amorphizat ion  of  materials  have  given 
rise  to  an  interest  in  the  point  defect  properties  of  mate¬ 
rials  that  can  be  amorphized,  such  as  NiTi  and  FeTi  (Sabo- 
chick,  1990,  in  press) . 

Recent  computer  simulation  studies  of  the  ordered  com¬ 
pound  copper  titanium  (CuTi)  metal  have  shown  that  the 
properties  of  point  defects  in  alloys  may  be  considerably 
different  than  those  of  their  pure  metal  counterparts 
(Shoemaker,  et  al.,  1990,  in  press).  Interstitials  in  this 
compound  were  found  to  have  complicated  configurations  con¬ 
taining  one  or  more  antisite  defects  (where  an  atom  is  on  a 
lattice  site,  but  is  of  the  wrong  type) .  The  complicated 
configuration  implies  the  interstitials  may  not  be  as 
mobile  as  they  are  in  pure  metals.  These  results  are  dif¬ 
ferent  from  what  is  observed  in  pure  metals  and  support 
recent  studies  of  the  amorphization  of  CuTi  (Shoemaker,  et 
al.,  1990,  in  press).  Computer  simulation  of  amorphization 
indicates  the  presence  of  point  defects  is  critical  in  FeTi 
lattices,  but  chemical  disorder  may  be  more  important  in 
NiTi  (Sabochick  and  Lam,  1990,  in  press).  A  comparison  of 
point  defect  properties  between  the  two  compounds  may 
reveal  the  reason  for  this  behavior. 
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Problem 


The  problem  is  lack  of  information  on  the  point  defect 
properties  of  the  intermetallic  compounds  nickel  titanium 
(NiTi)  and  iron  titanium  (FeTi) .  These  compounds  are 
important  because  they  can  be  amorphized  using  electron 
radiation.  In  addition,  because  NiTi  and  FeTi  have  the 
same  crystal  structure,  the  differences  in  their  defect 
properties  will  result  only  from  the  differences  between 
the  nickel  and  iron  elements.  Thus,  the  results  can  be 
directly  compared  for  understanding  the  effect  each  element 
has  on  point  defects. 

Scope 

The  problem  is  limited  to  single  atom  point  defects  in 
the  NiTi  and  FeTi  compound  lattices.  Thus  the  interaction 
of  the  defects  are  ignored,  and  the  theories  dependent  on 
small  concentrations  of  defects  can  be  used  to  explore  the 
defect  properties.  Specifically,  the  vacancies  and  inter¬ 
stitials  of  each  atom  type  in  each  compound  lattice  will  be 
explored  with  atomistic  simulation.  This  will  result  in  a 
total  of  eight  case  studies:  (1)  iron  vacancy,  (2)  iron 
interstitial,  (3)  titanium  vacancy  in  FeTi,  (4)  titanium 
interstitial  in  FeTi,  (5)  nickel  vacancy,  (6)  nickel  inter¬ 
stitial,  (7)  titanium  vacancy  in  NiTi,  and  (8)  titanium 
interstitial  in  NiTi. 
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Approach 


The  point  defect  properties  of  the  compounds  FeTi  and 
NiTi  are  studied  wich  computer  atomistic  simulation  tech¬ 
niques  with  a  modified  version  of  the  DYNAMO  code.  The 
general  approach  of  the  work  is  to: 

1.  Identify  the  stable  vacancy  configuration, 

2.  Calculate  the  migration  energy  of  the  vacancy, 

3.  Identify  the  stable  interstitial  configuration,  and 

4.  Attempt  to  calculate  the  migration  energy  of  the 
interstitial 

for  each  of  the  individual  compounds  NiTi  and  FeTi.  Inter¬ 
stitial  configuration  migration  of  similar  compounds  have 
proven  difficult  to  calculate  with  the  atomistic  simulation 
techniques  used  in  this  work.  The  stable  interstitial 
configuration  of  NiTi  and  FeTi  were  accompanied  by  antisite 
defects  (see  interstitial  results  section  V) .  Thus,  migra¬ 
tion  of  the  interstitial  was  not  pursued. 
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II .  Simulation  Method 


A  modified  version  of  the  DYNAMO  code  called  "RUNDYN" 
(see  Appendix  A  for  details)  was  used  in  the  atomistic  sim¬ 
ulation  required  for  the  calculations  in  this  work.  The 
simulation  techniques  included;  (1)  molecular  statics  for 
the  energy  migrations  and  minimizations,  (2)  molecular 
dynamics  with  simulated  annealing  for  vacancy  calculations, 
and  (3)  Monte  Carlo  atom  switching  for  the  interstitial 
calculations . 

In  molecular  dynamics  Newton's  equations  of  motion  for 
the  atoms  are  solved.  The  second  order  differential  equa¬ 
tions  are  solved  numerically  and  yield  the  position  and 
velocity  of  each  atom  at  discrete  times  for  the  duration  of 
the  simulation.  With  the  position  and  velocity  of  the 
atoms  the  lattice  properties  such  as  temperature,  pressure, 
and  internal  energy  are  calculated  (Sabochick,  1990,  in 
press) .  In  molecular  statics,  the  atomic  positions  are 
calculated  and  the  total  potential  energy  of  the  lattice  is 
minimized.  This  is  usually  done  with  mathematical  minimi¬ 
zation  techniques  such  as  steepest  descent  and  conjugate 
gradients.  In  the  modified  version  of  the  DYNAMO  code 
the  Fletcher-Powell  minimization  method  was  used  in  all  of 
the  calculations  of  this  work.  A  more  detailed  explana¬ 
tion,  with  regards  to  the  RUNDYN  code,  of  atomistic  simu- 


lation  is  discussed  by  M.  J.  Sabochick  in  his  report 
entitled  “Basic  Primer  on  Atomistic  Simulation." 

Additional  discussion  on  atomistic  simulation  is  found  in 
J.  M.  Haile's  report  entitled  "A  Primer  on  the  Computer 
Simulation  of  Atomic  Fluids  by  Molecular  Dynamics." 

All  of  the  calculations  had  cubic  volumes  with  periodic 
boundary  conditions  having  fixed  edge  lengths  and  constant 
pressure.  The  embedded  atom  method  (EAM)  was  used  with  the 
interatomic  potentials  derived  with  the  method  of  Oh  and 
Johnson  (Shoemaker,  et  al.,  1990,  in  press).  The  inter¬ 
atomic  potential  is  used  to  specify  the  energy  between 
atoms  as  a  function  of  distance.  These  potentials  were 
also  used  in  the  amorphization  studies  of  NiTi  and  FeTi  by 
M.  J.  Sabochick.  The  potentials  for  Ni,  Fe,  and  Ti  were 
fitted  to  the  binding  energies,  lattice  constants,  elastic 
constants  and  vacancy  formation  energy  of  the  pure  metals 
(Shoemaker,  et  al.,  1990,  in  press).  The  potentials  for 
FeTi  and  NiTi  were  then  determined  by  fitting  the  struc¬ 
ture,  lattice  constants  and  heat  of  formation  of  the  two 
compounds.  The  fitting  procedure  details  have  been 
described  by  Sabochick  and  Lam  in  other  works  (Sabochick 
and  Lam,  1990,  565) .  The  interatomic  potential  functions 
thus  become  the  foundation  for  the  point  defect  properties 
calculated  in  this  work. 
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The  potentials  used  give  lattice  parameters  for  FeTi  as 
2.976  +/-  0.005  Angstroms(A°  =  10'10  meters)  and  for  NiTi 
as  2.998  +/-  0.005  A°.  Both  of  these  values  are  slightly 
larger  than  experimental  values  given  by  Duwey  and  Taylor 
as  2.975  +/-  0.005  A°  for  FeTi  and  2.986  +/“  0.005  A°  for 
NiTi  (McQuillan,  1956,  221  and  239).  These  compounds  both 
have  body  B2-type  structures  described  as  two  interpene¬ 
trating  equal  sized  cubic  lattices  shown  in  Figure  1  (Gha- 
tak,  1972,  211).  Each  cubic  lattice  is  of  a  different  atom 
type  with  the  corner  atom  of  the  opposite  type  at  each 
cube's  center  position.  Thus,  in  this  simple  lattice 
arrangement  there  is  a  total  of  16  atoms  with  eight  atoms 
of  each  type. 


§ 


+  z 


NiTi  2.998  A 


Figure  1.  B2-type  Structure  of  FeTi  and  NiTi 
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Point  defects  (especially  interstitials)  have  many 
metastable  energy  configurations  that  require  a  technique 
to  determine  the  global  minimum  energy  (lowest  energy,  most 
stable)  configuration  (Shoemaker,  et  al.,  1990,  in  press). 
First,  a  lattice  structure  of  1024  atoms  was  created  with 
the  RUNDYN  code  for  both  NiTi  and  FeTi.  The  stable  and 
interstitial  global  minimum  energy  configurations  were  then 
determined  with  "simulated  annealing"  or  "Monte  Carlo  atom 
switching"  computer  techniques.  In  the  molecular  dynamic 
"simulated  annealing  technique,"  a  single  atom  was  either 
removed  (vacancy)  or  added  (interstitial)  at  or  near  the 
center  of  the  lattice  and  the  energy  of  the  lattice  was 
minimized.  These  vacancy  and  interstitial  lattices  were 
then  heated  to  800  degrees  Kelvin  (K)  and  minimized  to  the 
lowest  energy  configuration  using  the  Fletcher-Powell 
method  (Shoemaker,  et  al.,  1990,  in  press).  After  which, 
the  lattice  for  each  case  was  cooled  in  50  or  100  degree 
increments  to  about  room  temperature  at  300°  K.  The  lowest 
energy  configuration  (global  minimum)  of  these  calculations 
was  selected  for  each  case  from  each  group  of  calculations. 
Simulated  annealing  was  adequate  in  finding  the  lowest 
energy  configuration  of  vacancy  point  defects.  However, 
for  interstitial  point  defects  calculations  the  Monte  Carlo 
atom  switching  technique  proved  more  efficient. 


In  the  "Monte  Carlo  switching  technique,"  atoms  at  a 


constant  temperature  were  randomly  switched  from  one  type 
to  another  creating  aatisite  defects  (where  an  atom  is  on  a 
lattice  site,  but  is  of  the  wrong  type) .  The  total  number 
of  each  type  of  atom  is  conserved  in  the  switching  and  the 
lowest  energy  configuration  is  found  after  the  atoms  are 
switched.  This  technique  is  relatively  new  for  use  with 
the  RUNDYN  code  (Shoemaker,  et  al . ,  1990,  in  press)  and  was 
found  to  reduce  computational  time  and  to  be  more  efficient 
in  calculating  the  interstitial  configurations.  Both  the 
annealing  and  switching  techniques  were  used  in  searching 
for  the  interstitial  global  minimum  energy  configurations 
(most  stable) .  However,  the  final  configurations  were  usu¬ 
ally  found  with  a  little  additional  deductive  reasoning  to 
the  calculated  configurations. 

After  finding  the  lowest  energy  configuration  or  global 
minimum  the  results  of  the  vacancy  calculations  were 
explored  to  find  the  required  energy  for  individual  atoms 
to  migrate  between  lattice  sites.  Migration  energies  were 
calculated  by  minimizing  the  lattice  energy  as  a  jumping 
atom  (the  atom  moving  between  lattice  sites)  was  con¬ 
strained  at  various  distances  from  a  vacancy  in  the  lattice 
(see  Appendix  A  for  details) .  These  calculations  provided 
the  required  energy  for  atoms  to  move  within  the  lattice 
and  thereby  overcome  migration  energy  barriers. 


Error  assessment  is  not  generally  discussed  in  most 
atomistic  simulation  works.  Haile  makes  note  of  it  in  his 


report  concerning  the  numerical  solution  of  Newton's  equa¬ 
tions.  He  notes  that,  errors  are  "primarily  due  to  (a)  use 
of  a  finite-difference  algorithm  for  solving  differential 
equations  and  (b)  round-off  errors  that  occur  in  the  com¬ 
puter  hardware"  (Haile,  1980,  38).  In  all  of  the  calcula¬ 
tions  with  RUNDYN  seven  digits  after  the  decimal  point  are 
repeatable  if  the  tolerance  was  set  at  1.0  E-13  (see 
Appendix  B) .  The  results  are  thus  precise  but  without 
direct  experiment  results  may  lack  accuracy .  rnough  past 
work  has  been  done  with  the  DYNAMO  code  that  the  accuracy 
of  the  work  is  believed  to  be  within  a  few  tenths  (+/-  0.1 
eV)  of  an  eV  (Sabochick,  Dec  1990,  personal  discussion) . 

The  results  here  are  presented  with  3  or  4  significant  dig¬ 
its,  which  is  generally  adequate  for  comparing  results. 
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The  antisite  defect  (where  an  atom  is  on  a  lattice 
site,  but  is  of  the  wrong  type)  is  the  simplest  defect  of 
an  ordered  compound  such  as  NiTi  or  FeTi  (Figure  3) . 


> 


> 


+  z 


Ti  Antisite  Defect 


+z 


Ni  or  Fe  Antisite  Defect 


Figure  3.  Antisite  Defect  Lattice  Structure  for  NiTi  and 
FeTi 


The  formation  energy  of  the  antisite  defect  pair  is  easily 

>  calculated  with  the  RUNDYN  code  system  lattice  results. 

First,  a  single  atom  in  the  perfect  lattice  is  converted  to 
the  opposite  type  atom  (Figure  3)  and  a  minimization  calcu¬ 
li  lation  is  performed  on  the  resulting  lattice.  Second,  this 

is  done  for  each  atom  type  in  the  compound  lattice  and  the 
perfect  lattice  energy  (UQ)  is  subtracted  from  each  of  the 

>  antisite  defect  lattice  energies  (£/q‘  and  UTal ;  UFae  and  UTal)  . 
For  example,  in  the  NiTi  lattice  the  equations  for  the 
individual  antisite  defect  energies  (Zf£‘  and  ETa')  are: 

i 
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III.  Antisite  Defect 


The  perfect  lattice  structure  used  in  all  the  calcula¬ 
tions  was  a  cube  consisting  of  repeating  interpenetrating 
simple  cubic  structures  (Figure  1)  with  8  atoms  to  a  side 
having  side  dimensions  of  11.992  Angstroms (A°)  for  NiTi  and 
11.904  A°  for  FeTi  (Figure  2).  The  total  energy  (UQ)  of 
the  perfect  lattice  containing  1024  lattice  sites  (Ns)  was 
U  o‘  =  -5120.861  eV  for  the  NiTi  lattice  and  was  U  F0e  = 
-5008.453  eV  for  the  FeTi  lattice. 


NiTi=!  1  992(A)  or  FeTi= 1 1  904(A) - •» 


Figure  2.  Perfect  Lattice  Structure  used  for  NiTi  and  FeTi 
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(1) 


Ea=UNal~U  o 


E 


Ti 

a 


=  uTa'-u0 


(2) 


Finally,  the  Ni  antisite  pair  formation  energy  f/f  is 


calculated  by  adding  the  individual  antisite  defect  ener¬ 
gies  Eal  and  ETal. 


E  f1!'  =  E‘a+  ET‘ 


a 


(3) 


The  equations  (1),  (2),  and  (3)  are  identical  for  the  FeTi 

alloy  with  the  Ni  superscript  exchanged  with  a  Fe.  The 
antisite  pair  formation  energies  for  both  alloys  were 
calculated  from  the  code  results  given  in  Table  1  as  1.015 
eV  for  NiTi  and  0.726  eV  for  FeTi.  In  creating  the 
antisite  defects  the  energy  increase  was  larger  (less 
negative)  when  converting  Ti  atoms  to  either  nickel  (Ni)  or 
iron  (Fe)  than  when  converting  Ni  or  Fe  to  Ti  (Table  1) . 
This  implies  that  Ni  or  Fe  antisite  defects  should  be 
preferred  over  Ti  antisite  defects  in  the  lattice.  A  Ti 
vacancy  in  the  lattice  would  therefore  have  a  higher 
potential  of  being  filled  with  a  Ni  or  Fe  atom  and  thereby 
create  a  Ni  or  Fe  antisite  defect.  The  stable  vacancy 
configurations  given  in  the  next  section  confirm  this 
observation. 
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If  experimentally  neither  of  these  compounds  becomes 
chemically  disordered  when  heated,  it  would  imply  that  the 
antisite  defect  pair  formation  energy  is  significantly 
larger  than  the  Boltzmann's  constant  (k)  multiplied  by  the 
melting  temperature  (Tm)  of  the  given  alloy  (Shoemaker,  et 
al.,  1990,  in  press).  The  melting  temperatures  of  NiTi 
and  FeTi  are  1615°  K  and  1590°  K  respectively  (McQuillan, 
1956,  236  and  217).  Thus,  given  the  value  of  k  = 

8.6173  5E-5  eV  /  K  then  the  value  of  k*Tm  is  equal  to  0.139 
eV  and  0.137  eV  for  NiTi  and  FeTi  respectively.  These 
alloys  should  not  become  chemically  disordered  when  heated, 
because  the  antisite  formation  energies  are  significantly 
larger  (1.015  eV  for  NiTi  and  0.726  eV  for  FeTi)  then  their 
respective  k*Tm  values. 
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IV.  Vacancy  Defect  Configuration  and  Migration 


Ti  Vacancy  Defect  Properties 

A  single  Ti  atom  was  removed  from  the  center  of  the 
perfect  lattice  and  placed  at  the  end  of  the  restart  file 
for  both  NiTi  and  FeTi.  This  allowed  the  atoms  to  be 
renumbered  for  easily  tracking  the  movement  of  surrounding 
atoms.  This  same  Ti  atom  #1024  was  then  permanently 
removed  from  the  lattice  and  a  Ti  vacancy  was  created. 

Using  the  simulated  annealing  technique,  the  resulting 
minimum  lattice  energy  of  these  Ti  vacancy  lattices  was 
-5113.8369  eV  and  -5002.1175  eV  for  NiTi  and  FeTi  respec¬ 
tively.  The  stable  Ti  vacancy  configurations  had  higher 
energy  levels  than  the  perfect  lattice  by  7.0238  eV  for 
NiTi  and  6.3355  eV  for  FeTi  (Table  3).  Stability  of  these 
configurations  was  achieved  by  the  migration  of  an  adjacent 
Ni  or  Fe  atom  into  the  Ti  vacancy  (Part  A  in  Figure  4) . 

The  resulting  configuration  in  both  compounds  consisted  of 
a  vacant  Ni  or  Fe  site  and  an  adjacent  Ni  or  Fe  antisite 
defect  (Part  B  in  Figure  4) . 

The  migration  energy  barriers  to  the  migrating  Ni  or  Fe 
atoms  were  calculated  with  the  molecular  statics  technique 
described  in  the  simulation  section.  The  Ni  or  Fe  atom 
migration  energy  barriers  were  calculated  to  be  0.1344  eV 
for  NiTi  and  0.3407  eV  for  FeTi. 
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Figure  4.  Ti  Vacancy  Stable  Configuration  with  Ni  or  Fe 
Migration 


Each  of  these  migrations  is  shown  graphically  in  Figure  5 
for  equally  spaced  stepped  migration  distances.  See  Appen¬ 
dix  A  for  a  description  of  the  migration  method  and  Appen¬ 
dix  B  for  tabulated  migration  results.  In  the  final  stable 
configuration  the  adjacent  atoms  around  the  vacancy  were 
displaced  slightly  towards  the  vacancy.  This  relaxation 
around  the  vacancy  did  not  cause  large  displacements  of  the 
atoms  off  their  normal  lattice  sites  for  NiTi  nor  for  FeTi 
(see  Appendix  D  for  computer  output  plots) . 
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Figure  5.  Migration  Energy  Barriers  for  Ni  or  Fe  Atoms  to 
Ti  Vacancy 

The  Ni  or  Fe  antisite  defect  binding  energies  to  a  Ti 
vacancy  were  calculated  using  the  results  of  the  vacancy 
defect  calculations  given  in  Table  3  and  the  antisite 
defect  results  given  in  Table  1.  This  binding  energy  is 
the  energy  required  to  remove  the  Ni  or  Fe  antisite  defect 
from  the  lattice  and  replacing  the  resulting  vacancy  with  a 
Ti  atom.  This  was  done  by  adding  the  Ni  or  Fe  vacancy 
energies  to  the  Ni  or  Fe  antisite  defect  energies  for  their 
individual  alloy  types  (Table  2) .  The  Ti  vacancy  energies 
were  then  subtracted  from  the  above  results  for  their 


17 


respective  compounds.  This  resulted  in  the  binding  ener¬ 
gies  for  each  of  the  compounds  (Table  2) .  The  antisite 
defect  and  vacant  Ni  or  Fe  sites  were  found  to  be  weakly 
bound  with  a  binding  energy  of  0.094  eV  for  NiTi  and  0.055 
eV  for  FeTi. 


Table  2.  Ni  and  Fe  Antisite  Defect  Binding  Energies 

System 

Number  of  Atoms 

Energy  (eV) 

Ni  or  Fe 

Ti 

NiTi 

FeTi 

A.  Perfect  Crystal 

512 

512 

0.0 

0.0 

B.  Ti  Vacancy 

512 

511 

7.0238 

6.3355 

C.  Sum  =  A  +  B 

1024 

1023 

7.0238 

6.3355 

D.  Ni  or  Fe  Antisite 

513 

511 

0.7623 

0.4910 

E.  Ni  or  Fe  Vacancy 

511 

512 

6.3558 

5.8993 

F.  Sum  =  D  +  E 

1024 

1023 

7.1181 

6.3903 

Binding  Energy  =  F  -  C 

0.0943 

0.0548 

Ni  and  Fe  Vacancy  Defect  Properties 

The  Ni  or  Fe  vacancy  lattices  were  created  in  the  same 
way  as  the  initial  Ti  vacancy  lattices.  The  resulting 
relaxation  of  the  adjacent  atoms  surrounding  the  vacancy 
proved  to  be  the  lowest  energy  configuration,  as  no  read¬ 
justment  of  the  lattice  resulted  from  the  applied  simulated 
annealing  technique.  Thus,  the  stable  configuration  of  a 
vacancy  after  removal  of  a  Ni  atom  was  a  vacant  Ni  site  for 
both  compounds.  The  resulting  calculated  lattice  energy 
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values  are  given  in  Table  3 . 


Table  3.  Vacancy  Defect  Properties  of  NiTi  and  FeTi(£T  =  U*x-  Ua) 

NiTi  (y=Ni) 

FeTi  (y=Fe) 

System 

m 

tfXeV) 

fXeV) 

C?*eV) 

fXeV) 

Perfect  Crystal 

0 

-5120.861 

0.0 

-5008.453 

0.0 

Ni  or  Fe  Vacancy 

1 

-5114.505 

6.356 

-5002.554 

5.899 

Ti  Vacancy 

V  Ti 

-5113.837 

7.024 

-5002.118 

6.336 

The  migration  barrier  calculation  of  either  the  Ni  or  Fe 
vacancy  consists  of  a  six-jump  ring  sequence.  The  vacant 
Ni  or  Fe  site  is  rotated  through  its  nearest  neighbors 
(atoms  A  and  B  in  Figure  6)  to  the  next  Ni  or  Fe  lattice 
site.  The  migration  energy  barriers  calculated  through 
this  sequence  were  1.2718  eV  and  1.7382  eV  for  NiTi  and 
FeTi  respectively.  These  migration  values  are  relatively 
large  compared  to  typical  vacancy  migration  energies  in 
pure  metals.  This  can  be  attributed  to  the  intermediate 
configurations  containing  antisite  defects  and  or  vacant  Ti 
site  (Sabochick,  1990) . 
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Figure  6.  Ni  and  Fe  Vacancy  Six-Jump  Migration  Sequence 


The  migration  sequence  is  shown  graphically  in  Figure  7  and 
the  individual  migration  energies  are  given  in  Table  4. 

Both  the  Ni  and  Fe  migration  sequences  are  similar  except 
the  energy  barrier  occurs  in  the  second  and  fifth  steps  of 
the  sequence  for  the  FeTi  lattice  verses  the  third  and 
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Energy  (eV) 


forth  steps  for  the  NiTi  lattice  (Figure  7)  .  In  the  FeTi 
lattice  the  first  migration  sequence  also  results  in  an 
unstable  (no  metastable)  configuration. 


Equal  Spaced  Step  Migrations 


Figure  7.  Fe  and  Ni  Vacancy  Six-Jump  Migration  Graphic 
Results 

The  direct  migration  of  the  vacancy  to  its  second  nearest 
neighbor  location  was  also  calculated  for  comparison  to  the 
nearest  neighbor  (indirect)  migration.  These  direct  migra¬ 
tion  values  1.6594  eV  for  NiTi  and  2.3556  eV  for  FeTi  were 
significantly  higher  than  the  ring  sequence  migration 
sequence  for  both  compounds  (Figure  7  and  Table  4) . 
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Table  4.  Ni  and  Fe  Vacancy  Migration  Energies 


Migration  Path 


Ti  to  X  Vacancy 


X  to  Ti  Vacancy  (a) 


Ti  to  X  Vacancy  (b) 


Direct  Migration 


X  to  X  Vacancy 


(a)  with  1  adjacent  antisite  defect. 

(b)  with  2  adjacent  antisite  defects. 

(c)  unstable  configuration. 


X=Ni  (eV) 


X=Fe  (eV) 


1.206  (c) 
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V.  Vacancy  Defect  Formation  Energy 


Experiments  with  solids  in  which  point  defects  were 
observed  has  given  information  about  the  energies  of  defect 
formation,  volume  changes,  and  binding  energies  of  the 
defects.  From  this  information  the  theory  of  atomistic 
models  of  defects  has  developed.  The  theory  is  used  in 
atomistic  simulation  (1)  to  determine  the  configuration  of 
atoms  for  certain  (antisite,  vacancy,  or  interstitial) 
defects  in  a  given  material  and  (2)  to  calculate  the  values 
of  these  quantities  (Gruber,  1966,  21) . 

Experimental  quenching  techniques  have  been  used  in 
conjunction  with  direct  observation  by  ion  microscopy  to 
determine  the  energy  of  vacancy  defect  formation.  For 
example,  in  one  such  experiment  (Newkirk  and  Wernick,  1962, 
84)  a  platinum  tip  was  quenched  from  a  temperature  (7)  of 
1800°  K  and  five  vacancies  (N vac)  were  found  by  ion  micro¬ 
scopic  inspection.  The  total  number  of  atoms  (/V)  of  the 
inspected  material  was  8500  and  the  vacancy  concentration 
(C u  =  N„ac/N)  is  equal  to  5.9E-4  by  direct  counting.  Assum¬ 
ing  this  concentration  was  frozen  in  the  platinum  by  the 
quenching  and  assuming  the  vacancies  were  in  thermal 
equilibrium,  the  energy  of  formation  is  (E  f  =  kT  In  (C„)) 
equal  to  1.15  eV,  where  k  =  8.61735E-5  (qV  /  K)  is  the 
Boltzmann  constant.  This  technique  works  well  with  pure 
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metals  but  is  complicated  by  the  presence  of  other  atom 
types,  as  in  alloys.  Calculational  simulation  techniques 
can  however  be  used  determine  the  formation  energies  in 
alloys. 

In  pure  metals,  the  vacancy  is  of  only  one  type  of  atom 
without  the  possibility  of  antisite  defects.  The  formation 
energy  is  calculated  from  the  potential  energy  of  the  per¬ 
fect  lattice  and  the  energy  of  the  lattice  with  the  vacancy 
defect.  In  impure  metals  and  alloys  a  more  complex  variety 
of  point  defects  is  present,  these  include  antisite  defects 
and  vacancies  of  each  type  of  atom  in  the  alloy.  The  cal¬ 
culation  of  the  formation  energy  is  dependent  on  knowing 
the  potential  energies  01  each  of  these  point  defects.  One 
method  of  determining  the  vacancy  formation  energy  is  to 
calculate  the  concentration  of  the  defect  (C„)  as  a  func¬ 
tion  of  temperature  ( T )  where  the  defect  energy  of  forma¬ 
tion  (E f)  is  given  by 


Ef(T)=-k 


ri[ln(C,)] 

d[7-‘] 


(4) 


The  concentration  of  the  vacancy  defect  (C„)  as  a 


function  of  temperature  was  calculated  using  the  method  by 
Foiles  and  Daw  (Foiles  and  Daw,  1987,  12-13) .  In  this 
method,  the  concentrations  of  the  defects  are  assumed  small 
such  that  the  defects  can  be  treated  as  noninteracting. 
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For  the  NiTi  alloy,  there  is  one  Ni  lattice  site  for  each 
Ti  site.  Each  lattice  site  can  be  occupied  by  either  an 
atom  appropriate  to  the  sublattice,  by  a  vacancy,  or  by  an 
atom  of  the  opposite  type  (e.g.  antisite  defect) .  In  such 
a  system  for  the  NiTi  lattice  the  total  potential  internal 
energy  U  is 


U  =  if  N  E™  +  nTJ  El'  +  n%‘  E1?  +  nTa‘ETa‘)  <5) 


where  UQ  is  the  energy  of  the  perfect  lattice  with  Ns 
lattice  sites,  a"1,  nTvl ,  and  nTal  are  the  defect 

concentrations  of  Ni  vacancies,  Ti  vacancies,  Ni  antisite 
defects,  and  Ti  antisite  defects  respectively.  In  this 
system  the  total  concentration  of  the  vacancies  is  equal  to 
the  sum  of  the  individual  vacancies  of  each  type  of  atom 
(C „  =  n"1  +  nll)  .  The  corresponding  energies  E^\  ETV\  Ea\  and 
El1  are  the  differences  in  energy  between  the  system 
containing  a  single  defect  and  the  perfect  lattice  (e.g. 

E„l  =  U"1  -  U0)  .  These  energy  values  are  presented  for  the 
antisite  defects  in  Table  1  and  for  the  vacancy  defects  in 
Table  3.  The  configuration  entropy  S  of  the  system  is 


^s{^[s(2^:t)-s(2n^)]  +  i[s(2n:i)  +  s(2n^)] 


(6) 


where  s  is  the  ideal  entropy  function. 


s  (  x  )  =  —  [  a:  1  n  (  x  )  +  (  1  -  x  )  1  n  (  1  -  x  )  ]  ( 7 ) 
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The  total  number  of  Ni  atoms  (/vNt)  ,  is  given  by 


and  similarly  for  the  total  number  of  Ti  atoms  (A/Tl)  in  the 

configuration  with  the  Ni  and  Ti  superscripts  exchanged. 
The  antisite  concentration  was  shown  by  Foiles  and  Daw  to 
equal 

Nl  1  /  exp[-(f^  +  |iri-iiM‘)/A:T]  \  <9) 

2  \  1  +exp[-(£'i'l  +  |in-nw‘)/fcr]/ 

and  for  the  vacancy  concentration  to  equal 

Kt  if  exp[-(Ff  \  tl0> 

R "  2\  1  +exp[-(f^  +  (i‘v<)/fcT]/ 


where  n-Nl  and  [lTi  are  the  chemical  potentials  of  the 

respective  atom  types.  The  equations  for  nTa'  and  nJul  are 
identical  with  Ni  and  Ti  superscripts  exchanged  for  the 
respective  atoms.  The  chemical  potentials  can  be  elimi¬ 
nated  by  using  their  relationship  to  the  Gibb's  free  energy 
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U-kTS  =  m-  /V ,V1  +  m.  Ti  /V  Ti 


(ii) 


U  and  substituting  the  total  potential  energy  U  given  in 
equation  (5)  and  the  configuration  entropy  S  given  in 
equation  (6)  .  Requiring  the  total  number  of  atoms  at  a 
given  temperature  to  equal  the  stoichiometry  of  the 
compound 

1  NNi  (12) 

2  N  Ni  +  i V  Ti 

allows  the  eight  unknowns  N -Vl/  NTl,  n %l,  n"‘,  nTal ,  nTa' ,  i iNl, 

and  |iTl  to  be  solved  by  the  eight  equations  (8a, 8b), 

(9a, 9b),  (10a, 10b),  (11),  (12),  and  their  complements.  All 

of  these  equations  are  the  same  for  the  FeTi  alloy  with  the 
Ni  subscripts  exchanged  with  Fe  and  the  appropriate  Ti 
energy  values  used  for  FeTi.  These  equations  were 
implemented  into  the  TK  Solver  Plus  software  program  to 
facilitate  graphing  and  tabulating  the  results  (Appendix 
C) .  The  only  required  input  values  to  the  program  are  the 
five  total  system  energies  for  the  global  minimum  configu¬ 
rations  and  the  number  of  lattice  sites  for  each  type  of 
atom  in  the  compound.  These  energies  include  the  perfect 
crystal  U0  for  each  alloy,  the  antisite  defect  total 
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system  energies  for  NiTi  or  for  FeTi  uFae,UTa'  (Table 

1) ,  plus  the  vacancy  defect  total  system  energies  for  NiTi 
U*\UJul  or  for  FeTi  UF\Ul'  (Table  3). 

The  vacancy  concentrations  calculated  with  the  above 
system  of  equations  were  used  to  calculate  the  vacancy  for¬ 
mation  energy  of  each  compound.  This  was  done  by  calculat¬ 
ing  the  slope  of  the  log  of  C„  vs  1/7  shown  in  the  Figure  8 
Arrhenius  plot.  The  vacancy  formation  energy  is  a  function 
of  temperature,  however  at  lower  temperatures  (below  about 
1200°  K)  its  dependence  is  negligible.  The  curves  in  Fig¬ 
ure  8  are  essentially  straight  indicating  that  a  single 
formation  energy  value  can  be  used  over  the  entire 
temperature  range.  The  effective  vacancy  formation  ener¬ 
gies  were  calculated  to  be  1.48  eV  for  NiTi  and  1.07  eV  for 
FeTi.  These  values  compare  reasonably  well  with  the  values 
of  CuTi  (1.09  eV)  and  CuTi2  (r'.3C  c7)  ,  -.."hid.  were  calcu¬ 
lated  with  the  same  method  from  data  given  by  Shoemaker's 
DYNAMO  calculations  (Shoemaker,  et  al.,  1990,  in  press). 

No  experimental  values  for  the  compounds  were  found,  how¬ 
ever  these  formation  energies  are  generally  lower  then 
those  of  the  individual  pure  metals.  For  example,  the 
experimental  formation  energies  of  the  metals  Ni  (1.7  +/- 
0.1  eV,  Agullo-Lopez  et  al.,  1988,  197)  and  Ti  (1.55  +/- 
0.1  eV,  Shoemaker,  et  al.,  1990,  in  press)  are  larger  than 
those  of  the  NiTi  compound  (1.48  +/-  eV) . 
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Vacancy  Defect  Concentration  (Cv) 


Inverse  Temperature  1000/T  (K~-1) 


Figure  8.  Arhennius  Plot  of  Vacancy  Concentration  vs 
perature 


VI . 


The  interstitial  calculations  had  the  same  perfect  lat¬ 
tice  structure  as  the  vacancy  calculations.  A  single  tita¬ 
nium  (Ti)  atom  was  added  near  the  center  of  the  lattice  by 
placing  it  at  the  end  of  the  restart  file.  This  again 
allowed  the  atoms  movement  to  be  easily  tracked  by  their 
atom  numbers  and  avoided  graphical  viewing  problems  asso¬ 
ciated  with  the  edge  of  the  lattice.  The  resulting  mini¬ 
mized  lattice  became  the  basis  for  all  the  other 
calculations . 

The  annealing  simulation  technique  proved  to  be  an 
ineffective  method  for  finding  the  global  minimum  configu¬ 
ration.  It  appears  the  annealing  technique  did  not  allow 
the  antisite  defects  to  be  formed  easily.  The  Monte  Carlo 
switching  technique  easily  formed  antisite  defects  and 
proved  to  be  the  most  useful  in  determining  the  intersti¬ 
tial  configuration.  However,  in  all  cases  the  final  con¬ 
figuration  was  determined  with  a  little  added  inductive 
reasoning  based  on  the  Monte  Carlo  calculation  results. 

The  Monte  Carlo  technique  is  dependent  on  finding  a  criti¬ 
cal  temperature  to  allow  switching  to  occur  in  the  lattice 
without  overwhelming  the  process  with  to  many  antisite 
defect  switches.  For  FeTi  no  switches  occurred  at  520 
degrees  K  and  too  many  occurred  at  550  degrees  K.  Thus  the 
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critical  temperature  used  for  FeTi  was  530  K  and  for  NiTi 
was  720  K.  Once  the  critical  temperature  was  established 
the  calculation  could  be  repeated  with  different  random 
seed  numbers,  which  should  generate  similar  results.  All 
of  the  interstitial  configurations  were  generated  in  this 
same  fashion. 

Ti  Interstitial  Defect  Properties 

The  final  stable  configuration  consisted  of  a  Ni-Ni 
or  Fe-Fe  dumbbell  centered  on  a  Ti  site  with  2  adjacent  Ti 
antisite  defects  (Ti  on  Ni  or  Fe  site)  at  opposite  corners 
of  the  Ni  or  Fe  cube  (Figure  9  and  10) .  These  configura¬ 
tions  were  the  same  for  NiTi  and  FeTi  except  for  the  direc¬ 
tion  of  the  dumbbell.  The  Ni-Ni  dumbbell 

(split-interstitial)  was  oriented  in  the  <111>  unit  vector 
direction  with  relatively  large  displacements  (off  the  lat¬ 
tice  site)  of  the  atoms  along  the  <111>  direction.  Output 
plots  of  these  configurations  are  given  in  Appendix  D. 
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Figure  9.  Ti  Interstitial  Stable  Configuration  for  NiTi 

+  z 


Figure  10.  Ti  Interstitial  Stable  Configuration  for  FeTi 

No  quantitative  measurements  of  the  displacements  were  made 
and  the  subjective  terms  large,  small,  or  slight  are  rela¬ 
tive  to  each  other.  RUNDYN  code  output  plots  of  all  of  the 
configurations  are  given  in  Appendix  D  for  comparison  of 
displacements.  The  Fe-Fe  dumbbell  was  oriented  in  the 
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<011>  direction  (Figure  10)  with  significant  but  smaller 
displacement  of  the  neighboring  atoms  in  the  <011?  direc¬ 
tion.  The  values  of  the  stable  configuration  energies  and 
Frenkel  pair  formation  energies  are  given  in  Table  4. 


Table  5.  Interstitial  Defect  Properties  (F(  =  Uyx-U0) 

NiTi  (y=Ni) 

FeTi  (y=Fe) 

System 

m 

u%eV) 

F^eV) 

v’JLeV) 

FXeV) 

Perfect  Crystal 

0 

-5120.861 

0.0 

-5008.453 

0.0 

Ni/Fe  Interstitial 

i  y 

-5123.256 

-2.395 

-5010.011 

-1.558 

Ti  Interstitial 

i  Ti 

-5123.463 

-2.602 

-5010.398 

-1.945 

Ni/Fe  Frenkel  Pair 

am 

3.961 

4.341 

Ti  Frenkel  Pair 

f  Ti 

4.422 

4.391 
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Ni  and  Fe  Interstitial  Defect  Properties 


The  Ni  and  Fe  interstitial  configurations  were  calcu¬ 
lated  in  the  same  manor  as  the  Ti  interstitial.  The 
resulting  global  minimum  energy  configurations  (Figure  11) 
for  both  NiTi  and  FeTi  consisted  of  a  Ni-Ni  or  Fe-Fe  dumb¬ 
bell  oriented  in  the  <111>  unit  vector  direction  with  an 
adjacent  single  Ti  antisite  defect  (Ti  on  Ni  or  Fe  site) . 
The  stable  defect  configuration  energies  and  calculated 
point  defect  energies  are  given  in  Table  4  for  all  the 
interstitial  defects. 
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Figure  11.  Ni  or  Fe  Interstitial  Configuration  for  NiTi 
and  FeTi 
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VII. 


Discussion 


The  stable  antisite  defect  configurations  (Table  6) 
were  calculated  and  used  to  determine  the  formation  ener¬ 
gies  and  binding  energies  (Table  7)  for  both  NiTi  and  FeTi. 
The  resulting  NiTi  antisite  pair  formation  energy  is  over 
twice  as  that  of  the  energy  for  either  CuTi  or  CuTi2- 
Sabochick  and  Lam  explored  the  key  factors  for  inducing  the 
crystalline-to-amorphous  (C-A)  transition  in  NiTi.  They 
concluded  that  point  defects  are  necessary  for  amorphiza- 
tion  to  occur,  but  that  chemical  disorder  was  sufficient  to 
store  the  energy  necessary  to  make  the  transition 
(Sabochick  and  Lam,  1990,  in  press) .  The  large  energy  nec¬ 
essary  for  NiTi  antisite  pair  formation  explains  why  chemi¬ 
cal  disordering  can  store  large  amounts  of  energy  in  NiTi. 
With  the  large  FeTi  antisite  pair  energy,  one  would  expect 
similar  amorphization  behavior  in  FeTi  as  is  found  in  the 
NiTi  alloy. 

The  stable  vacancy  configurations  of  both  NiTi  and  FeTi 
preferred  Ni  or  Fe  sites  for  vacancies  in  both  the  Ti 
vacancy  and  Ni/Fe  vacancy  calculations.  The  Ti  vacancy  was 
filled  by  a  migrating  Ni  or  Fe  and  resulted  in  a  Ni  or  Fe 
antisite  defect  with  a  neighboring  vacant  Ni  or  Fe  site. 
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This  same  site  preference  was  noted  to  occur  in  both  the 
CuTi  and  CuTi2  compounds  (Shoemaker,  et  al.,  1990,  in 
press)  . 


Table  6.  Stable  Defect  Configuration  Energies  (E  =  U 

-  U0) 

System 

Number  of 

NiTi  (eV) 

FeTi  (eV) 

Atoms 

Ti  Antisite  (a) 

■HUH! 

msm 

MEM 

Ni  or  Fe  Antisite  (b) 

1 

■9 

H 

Ti  Vacancy 

MEM 

7.024 

6.336 

Ni  or  Fe  Vacancy 

6.356 

5.900 

Ti  Interstitial  (c  &  d) 

msm 

-2.602 

MEM 

Ni  or  Fe  Interstitial  (e) 

KB 

-2.395 

WBEM 

(a)  Ti  on  Ni  or  Fe  Lattice  Site. 

(b)  Ni  or  Fe  on  Ti  Lattice  Site. 

(c)  NiTi  <11 1>  split-interstitial  with  two  Ti  antisite  defects. 

(d)  FeTi  <01 1>  split-interstitial  with  two  Ti  antisite  defects. 

(e)  Both  <U1>  split-interstitials  with  one  Ti  antisite  defect. 

It  is  well  accepted  that  vacancy  migration  is  the 
dominant  mechanism  for  self-diffusion  in  metals  (Agullo- 
Lopez  and  others,  3988,  196).  The  vacancy  migration 
energies  (Table  7)  for  NiTi  (1.27  eV)  and  FeTi  (1.74  eV) 
are  similar  to  those  for  CuTi  (0.19/1.32  eV)  and  CuTi2 
(0.92).  The  higher  FeTi  energy  indicates  diffusion  is  more 
difficult  in  this  alloy  than  either  NiTi,  CuTi,  or  CuTi2« 
The  higher  vacancy  migration  energies  may  be  explained  by 
(1)  the  antisite  defects  formed  in  the  migration  sequence 
and  by  (2)  the  large  antisite  pair  energies  of  NiTi  and 
FeTi.  The  vacancy  formation  energy  for  FeTi  (1.07  eV)  is, 
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however,  significantly  lower  than  that  of  NiTi  (1.48  eV) . 
Thus,  although  vacancies  are  more  easily  formed  in  FeTi 


their  migration  is  rore  difficult.  Although  no  experimen¬ 
tal  vacancy  formation  energies  were  found  for  the  alloys, 
these  calculated  values  are  lower  than  the  experimental 
values  for  the  pure  metals  of  the  compounds.  For  example 
for  NiTi  the  pure  metal  values  are  1.7  (+/-  0.1)  eV  for  Ni 
(Agullo-Lopez  and  others,  1988,  197)  and  1.55  eV  for  Ti 
(Shoemaker,  et  al.,  1990,  in  press).  This  indicates  the 
vacancy  formation  is  easier  in  both  alloys  NiTi  and  FeTi 
than  in  their  pure  metals.  The  same  results  was  noted  by 
Shoemaker  for  CuTi  and  CuTi.2- 


Table  7.  Calculated  Point  Defect  Energies  (in  eV) 

System 

NiTi 

FeTi 

CuTi  (a) 

CuTi2  (a) 

Antisite  Pair  Formation 

1.015 

0.726 

0.385 

0.460 

Antisite  Binding 

0.094 

0.055 

Vacancy  Migration 

1.272 

1.738 

0.19/1.32 

0.92 

Vacancy  Formation 

1.48 

1.07 

1.09 

0.90 

Frenkel  Pair  Ti  Formation 

4.422 

4.391 

Frenkel  Pair  Ni,  Fe,  or  Cu 

3.961 

4.341 

2.769 

2.776 

Interstitial  Formation  (b) 

2.5 

3.3 

1.7 

1.9 

(a)  Shoemaker,  et  al.,  1990,  in  press. 

(b)  Estimated  by  subtracting  the  Vacancy  Formation  from  the  Frenkel  Pair 
Energy. 

The  interstitial  configurations  of  both  NiTi  and  FeTi 
were  identical  for  Ni  or  Fe  interstitials  and  almost 
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identical  for  Ti  interstitials.  These  configurations  all 
involved  split-interstitials  with  neighboring  antisite 
defects.  The  antisite  defects  increase  the  migration 
energies  of  these  alloys,  because  of  the  large  number  of 
atoms  involved  in  the  migration  (the  split-interstitials 
and  the  antisite  defects) .  The  associated  antisite  defects 
with  the  interstitial  configurations  also  make  migration 
using  molecular  statics  difficult.  The  dumbbell  arrange¬ 
ment  as  well  as  the  associated  antisite  defects  must  be 
migrated  together  or  in  some  reasonable  order  to  achieve 
migration  of  the  defect  as  a  unit.  Since  the  method  of 
this  migration  was  not  obvious,  migration  of  the  intersti¬ 
tial  defect  was  not  done  for  either  NiTi  or  FeTi. 

An  estimate  of  the  interstitial  formation  energy  can  be 
made  from  the  Frenkel  pair  formation  energies  and  vacancy 
formation  energies  for  each  alloy.  The  Frenkel  (a  vacancy 
and  interstitial  defect  combination)  pair  formation  ener¬ 
gies  for  FeTi  are  4.34  eV  for  the  Fe  pair  and  4.39  eV  for 
the  Ti  pair  (Table  7) .  The  closeness  of  these  values  helps 
explain  the  similarity  of  the  vacancy  and  interstitial  con¬ 
figurations.  Subtracting  the  vacancy  formation  energy  from 
the  Frenkel  pair  gives  an  estimate  of  the  interstitial 
formation  energy.  The  vacancy  formation  energy  from  Table 
7  for  FeTi  is  1.07  eV  and  the  interstitial  formation  energy 
is  estimated  to  be  between  3.2  eV  and  3.3  eV.  For  NiTi 
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the  Frenkel  pair  energies  are  3.96  eV  for  the  Ni  pair  and 
4.42  eV  for  the  Ti  pair.  The  vacancy  formation  energy  is 
1.40  and  the  estimated  interstitial  formation  energy  is 

between  2.5  eV  and  2.9  eV.  This  implies  that  interstitial 
defect  formation  in  the  FeTi  lattice  is  more  difficult 
(higher  energy)  than  in  the  NiTi  lattice  with  a  difference 
in  energy  formation  of  0.3  to  0.8  eV. 

The  interstitial  configurations  with  a  <111>  oriented 
Ni-Ni  or  Fe-Fe  dumbbell  for  NiTi  and  FeTi  are  similar  to 
those  found  for  CuTi  and  CuTi2  (Shoemaker,  et  al.,  1990,  in 
press) .  Both  the  CuTi  and  CuTi2  alloys  have  a  Cu-Cu  dumb¬ 
bell  (split-interstitials)  centered  on  a  Ti  lattice  site  in 
the  <111>  unit  vector  direction  with  associated  neighboring 
antisite  defects.  The  exception  to  this  configuration  is 
the  Fe  interstitial  configuration,  because  it  is  oriented 
in  the  <011>  direction  verses  <111>  for  the  other  intersti¬ 
tials. 
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VIII.  Conclusion 


Atomistic  simulation  with  a  modified  version  of  the 
DYNAMO  code  was  used  to  calculate  the  point  defect  proper¬ 
ties  of  the  intermetallic  compounds  NiTi  and  FeTi.  The 
calculated  energies  are  believed  to  be  experimentally 
accurate  to  within  +/-  0.1  eV. 

The  antisite  pair  formation  energies  for  both  alloys 
were  calculated  to  be 

1.015  eV  for  NiTi  and  0.726  eV  for  FeTi.  These  values  are 
more  than  twice  the  calculated  formation  energies  of  CuTi 
(0.385  eV)  and  CuTi2  (0.460  eV) .  The  high  formation  ener¬ 
gies  in  NiTi  explains  why  large  amounts  of  energy  are 
required  for  chemical  disordering  in  the  compound.  The 
antisite  defect  and  vacant  Ni  or  Fe  lattice  sites  were 
found  to  be  weakly  bound  with  a  binding  energy  of  0.094  eV 
for  NiTi  and  0.055  eV  for  FeTi. 

The  stable  vacancy  defect  configurations  in  both  NiTi 
and  FeTi  were  identical.  Vacant  Ni  or  Fe  lattice  sites 
were  stable  and  preferred  in  the  alloys  with  energy  values 
of  6.356  eV  for  NiTi  and  5.899  eV  for  FeTi.  Ni  and  Fe 
vacancy  migrations  follow  a  six-jump  ring  sequence  with 
energy  barrier  values  calculated  as  1.272  eV  for  NiTi  and 
1.737  eV  for  FeTi.  Vacant  Ti  sites  resulted  in  neighboring 
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antisite  defects  from  migrating  Ni  or  Fe  atoms  with  energy 
values  of  7.024  eV  for  NiTi  and  6.336  eV  for  FeTi.  The 
migrating  Ni  ^>r  Fe  atoms  to  the  Ti  vacancies  had  migration 
barriers  values  of  0.134  eV  for  NiTi  and  0.341  eV  for  FeTi. 
The  vacancy  defect  formation  energies  were  calculated  from 
vacancy  concentration  values  to  be  1.48  eV  for  NiTi  and 
1.07  eV  for  FeTi. 

The  most  stable  interstitial  defect  configurations  con¬ 
sisted  of  Ni-Ni  or  Fe-Fe  split-interstitial  dumbbells  cen¬ 
tered  on  a  Ti  site  with  one  or  two  adjacent  Ti  antisite 
defects.  The  only  difference  between  the  two  alloys  was  in 
the  direction  of  the  Ti  interstitial  dumbbell,  whose  Fe-Fe 
orientation  was  in  the  <011>  direction.  Ti  interstitials 
resulted  in  two  adjacent  Ti  antisite  defects  with  a  <011> 
direction  Fe-Fe  dumbbell  or  a  <111>  direction  Ni-Ni  dumb¬ 
bell.  The  Ti  interstitial  formation  energies  were  -2.395 
eV  for  NiTi  and  -1.558  eV  for  FeTi.  Ni  or  Fe  interstitials 
both  resulted  in  a  single  adjacent  Ti  antisite  defect  with 
Ni-Ni  or  Fe-Fe  dumbbells  oriented  in  the  <111>  direction. 
The  Ni  or  Fe  interstitial  energies  were  -2.602  eV  for  NiTi 
and  -1.945  eV  for  FeTi.  Frenkel  pair  energies  for  NiTi 
were  3.961  and  4.422  eV  for  the  Ni  and  Ti  pair  respec¬ 
tively.  For  FeTi  the  Frenkel  pair  formation  energies  were 
4.341  and  4.391  eV  for  Fe  and  Ti  pairs  respectively. 
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APPENDIX  A: 


The  following  is  a  sample  copy  of  an  input  file  used  in 
a  molecular  statics  minimization  calculation.  A  descrip¬ 
tion  of  the  RUNDYN  code  is  also  included  for  understanding 
of  the  code's  basics  (Sabochick,  1990,  unpublished).  Since 
the  code  is  constantly  being  added  to,  this  information  is 
only  a  description  of  the  RUNDYN  version  used  for  this  work 
at  this  time. 

THE  CODE 

The  source  file  contains  a  computer  code  for  doing 
molecular  dynamics  and  molecular  statics  calculations  using 
embedded-atom  type  potentials.  The  code  is  a  modified  ver¬ 
sion  of  the  DYNAMO  code,  distributed  by  the  Sandia  group, 
with  the  following  modifications: 

1)  The  potential  files  (containing  tables  of  the  values 
for  the  components  of  the  EAM  potentials)  have  a  different 
format. 

2)  The  potentials  are  fitted  within  the  code  with  cubic 
splines  instead  of  Lagrange  polynomials. 

3)  The  code  can  do  molecular  statics  calculations  using  a 
modified  Fletcher-Powell  minimization  technique,  which  is 
faster  than  conjugate-gradients  on  large  problems. 
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4)  The  code  can  now  calculate  migration  energy  barriers  in 
a  scraight-f orward  manor. 

5)  Other  special  options  have  been  added  such  as;  atom 
switching,  Frenkel  pair  generation,  and  Monte  Carlo  anti¬ 
site  defect  switching. 


MIGRATION  ENERGY  TECHNIQUE 

Migration  energies  are  calculated  by  minimizing  the 
system  energy  using  the  Fletcher-Powell  method  with  a  jump¬ 
ing  atom  constrained  a  various  points  along  the  migration 
path  (Shoemaker,  et  al.,  1990,  in  press).  A  reaction 
coordinate  u  is  defined  so  that  the  points  r(r|)  along  the 
migration  path  from  r,  to  r2  are  specified  by 

r(ri)=ri  +  r)[r2-r1]  (13) 


The  lattice  system  is  then  minimized  with  the  jumping  atom 
confined  to  a  plane  perpendicular  to  r2-r,  at  the  point 
r(r|).  The  migration  barrier  energy  E(t])  is  then  determined 
by  repeating  the  calculation  for  various  values  of  r|,  and 
the  migration  energy  Em  is  the  maximum  of  E(r|)  or  rj  in  the 
range  -1  to  1  (i.e.  -1.0,  -0.75,  -0.5,  -0.25,  0.0,  0.25, 
0.5,  0.75,  and  1.0). 
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RUNNING  THE  CODE 


The  source  code  needs  to  be  compiled  and  linked  to  a 
input  file  (in  this  case  named  feti.inp).  A  softlink  was 
made  to  the  RUNDYN  source  code  in  M.  J.  Sabochick's  file 
directory.  The  input  files  could  then  be  run  by  typing 
RUNDYN  and  the  filename  of  the  input  file  (i.e.  feti.inp). 
The  output  file  name  (which  is  specified  in  the  input  file) 
is  feti.prn  for  this  example. 

INPUT  FILE  DESCRIPTION 


The  following  is  a  breakdown  of  the  input  file  "fe¬ 
ti.inp".  The  lines  starting  with  *****  and  in  bold  print 
are  the  input  file  lines  and  the  others  are  comments.  At 
the  beginning  of  each  section  is  a  header  card  of  the  form 
$xyz  which  helps  the  user  delimit  the  sections.  In  the 
newest  code  version  "RUNDYN2"  the  header  cards  and  their 
contents  can  be  omitted  when  the  section  is  not  used.  In 
the  previous  version  these  header  cards  were  checked  by  the 
code  (routine  chkname ( ) )  to  make  sure  it  was  in  the  right 
place.  Many  of  the  data  values  in  the  input  file  are  fol¬ 
lowed  by  a  description  (usually  the  name  of  the  variable, 
which  is  set) .  These  names  are  there  just  to  help  the  user 
and  are  optional.  All  input  lines  start  in  column  1. 
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PRINT  CONTROL  OPTION: 

*****  $prntcard 

Obsure  options  (not  used  in  the  calculations) 

*****  0  ipinter 

*****  0  ipatoms 

*****  -1  ipitera 

Name  of  output  file. 

*****  feti.prn 

Name  of  restart  save  file  (can  be  used  to  start  another 
run)  . 

*****  feti.res 

Name  of  the  trajectory  save  file. 

*****  none 

Name  of  the  energy  information  save  file. 

*****  none 


HEADER  INFORMATION: 

*****  $headcard 

Specific  Run  name  for  this  calculation. 
*****  Setup  of  FeTi  lattice 


ENERGY  POTENTIAL  FOR  RUN; 

The  potentials  are  referred  to  by  numbers  within  the  code; 
these  numbers  correspond  to  atom  types. (i.e.  1  for  Ni  or  Fe 
and  2  for  Ti) 

*****  $funccard 

Number  of  different  potentials  to  be  read  in. 

*****  2  ntypes 

File  Location  &/Name(s)  of  the  potential  files  (one  on  each 
line. ) 

*****  /enfacl/enp/msabochi/Johnson/Pot/fe03 .pot 

*****  /enfacl/enp/msabochi/Johnson/Pot/ti20 .pot 

*****  /enfacl/enp/msabochi/Johnson/Pot/feti_05 .pot 


The  code  can  start  from  a  restart  file  written  by  a  pre¬ 
vious  run,  or  generate  a  new  lattice. 

*****  $initcard 

Logical - whether  or  not  to  generate  a  lattice. 

*****  .true.  genlat 

Logical - whether  or  not  to  initialize  velocities. 

*****  .true.  genvel 

Logical - whether  or  not  to  sort  atoms  according  to  dis¬ 

tance  from  the  origin. 

*****  .false.  sort 

Logical - name  of  start  file  where  atom  posi^ions/veloci- 

ties  are. 

*****  none 
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used  if  genlat  equals 


true)  . 

*****  $latcard 

Size  of  system  in  unit  cells  (x,y,z). 

*****  888  ncells 

Unit  cell  dimensions.  Diagonal  matrix  is  a  cubic  or  rect¬ 
angular  cell. 

*****  2.976  0.000  0.000  basis  (fitted  a, 

a=2 .976) 

*****  0.000  2.976  0.000 

*****  0.000  0.000  2.976 

Lower  system  boundaries.  (x,y,z). 

*****  0.0  0.0  0.0  perlb 

Upper  system  boundaries.  If  values  are  <  0,  they  are  auto¬ 
matically  calculated  from  the  information  above. 

*****  -1.0  -1.0  -1.0  perub 

Number  of  atoms  in  a  unit  cell. 

*****2  nc 

Type  &  positions  of  atoms  in  unit  cell,  in  unit  cell  coor¬ 
dinates  (x,y , z) . 

*****  10. 00. 00. 0  NiorFe  type  1 

*****  2  0.5  0.5  0.5  Ti  type  2 


VELOCITY  INITIALIZATION  INFORMATION:  (only  used  if  genvel 
equal  true) . 

*****  $velcard 

Temperature  (Kelvin)  to  which  the  velocities  are  ini¬ 
tialized. 

*****  160.  temp 


SYSTEM  BOUNDARY  INFORMATION; 

*****  $bndcard 

Boundary  type.  (l=fixed  boundaries,  2=flexible,  rectangu¬ 
lar  boundaries) . 

*****  2  ibdtype 

The  following  six  variables  are  only  used  when  ibdtype 

equal  2: 

External  applied  pressure  (in  bars). 

*****  0.0  dpress 

The  following  five  variables  are  only  used  in  dynamic  cal¬ 
culations  ($intcard) :  Boundary  mass. 

*****  1.  i.  i.  bndmas 

Drag  (Automatically  set  if  negative) . 

*****  -1.  -i.  -i.  bnddrg 

Temperature  setting  option. 

*****  1 

Desired  Temperature  (Kelvin) . 

*****  160. 

Relaxation  Time. 

*****  0.100 
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ifxtmp 

destmp 

tmptim 


whether  or  not  to  use 


Option  to  use  Images. 
*****  .false, 

images 


NEIGHBOR  OPTIONS: 

*****  $neicard 

Neighbor  list  used?  (l=no,  2=yes) . 

*****  2  nmeth 

Cutoff  radius  (Angstroms)  for  neighbor  list  (only  used  if 
nmeth  equal=2) . 

*****  l.o  dradn 


DEFECT  INFORMATION:  (move 

cancies) . 

*****  $defcard 

Defect  type.  (-l=vacancy, 
99=not  used) . 

*****  99 

Atom  position  (x,y,z). 
*****  0.0  0.0  0.0 
Atom  velocity  (x,y,z). 
*****  0.0  0.0  0.0 


atoms,  create  interstitials/va- 

0=end  of  list,  l=interstitial , 
ktype 
pos 
vel 


Atom  number  to  delete  for  vacancy. 
*****  0  num 


TRAJECTORY  INFORMATION: 

*****  $trajsave 

*****  -1  save  interval 

Option  to  save  position  or  position  &  velocity. 

*****  1  what  is  saved  (l=pos 

only,  2=pos  and  vel) 


FORCED  "STATIC  MIGRATION"  MODIFICATIONS:  (added  to  deter¬ 
mine  migration  energy) . 

*****  $modcard 

Option  (0=0ff,  l=On) 

*****  o 

Number  of  Neighboring  Gate  atoms. 

*****  6 

First  #  is  Migration  atom  the  remaining  are  a  list  of  gate 
atom  numbers. 

*****  440  568  456  442  457  569  583 

Percent  Distance  r|(+or-)  from  Midpoint  r(r|)  (0.0),  Distan¬ 
ce  (A°),  Constant. 

*****  -0.25  1.2582  1.0e+4 

Unit  Vector  Direction  (x'-x,y '-y, z '-z)  from  Position 
(x,y,z)  to  (x ' , y ' , z ' ) • 

*****  1.0  1.0  1.0 
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MONTE  CARLO  ATOM  SWITCHING  INFORMATION:  (used  in  dynamic 
calculations) 


*  *  4t 

$swicard 

***** 

—  n 

kbsvit 

***** 

-2 

Antisite 

Defect  Energy 

(eV) 

Fe 

Atom  on  Ti 

Site. 

***** 

0.4910 

Antisite 

Defect  Energy 

(eV) 

Ti 

Atom  on  Fe 

Site. 

*****  0.2352 


FRENKEL  PAIR  CREATION  INFORMATION: 
*****  $frencard 

Option  (-l=Of f ) 

*****  — 1 

*****  0.7 

*****  i 

*****  2. 

*****  5.0 

*****  10.0 


SPECIFICATION  OF  RUN  TYPE:  (dynamic,  static) 

*****  $intcard 

Type  of  Integrator/Minimizer . 

(  1  =  Dynamic  Gear  predictor-corrector) 

(-1  =  Static  Conjugate-gradients  Minimizer) 

(-2  =  Static  Fletcher-Powell  Minimizer) 

(1000  =  Dynamic  Monte  Carlo  Switching  &  Pred/Corr) 


*****  -2 

inte 

Equilibrium  Steps 
*****  300 

(Dynamic 

runs) . 

equilibrium  steps 

Property  Steps 
*****  5000 

(Dynamic 

runs) . 

property  steps 

Output  Interval 
*****  10 

(Dynamic 

runs) . 

print  interval 

Accuracy  Option 
*****  o 

(Dynamic 

runs) . 

iaccur 

Time  Step 
*****  0.001 

(Dynamic 

runs) . 

dt 

Convergence  tolerance  (Static  Minimization  runs) . 
*****  l.e-13  tol 

Maximum  number  of  force  evaluations. 

*****  100000  nfmax 
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APPENDIX  B: 


Calculation  Spreadsheet  Results 


The  following  tables  are  a  copy  of  a  Lotus  123  spread¬ 
sheet  file  (MIGRATE. WK1)  used  to  record  and  manipulate  the 
calculated  lattice  configuration  energies.  This  data  was 
used  to  create  the  various  graphs  and  tables  given  in  the 
text.  The  system  energies  are  given  untruncated  for  com¬ 
pleteness  and  although  seven  digits  are  given,  only  3  are 
considered  significant  (Sabochick,  1990,  personal 
conversion) -  In  repeat  calculations  seven  significant  dig 
its  could  be  duplicated,  if  the  tolerance  was  input  as 
1.0E-13  (see  Appendix  A). 
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Table  8 


FeTi  "Ti"  Vacancy  MIGRATION  (Min 


-5002.1175360) 


Table  9 


%Distance 

Energy  (eV) 

Difference 

0.0000000 

-5001.0895532 

1.0279828 

0.1250000 

-5000.9724442 

1.1450918 

0.2500000 

-5000.7488333 

1.3687027 

0.3750000 

-5000.7767544 

1.3407816 

0.5000000 

-5001.0400101 

1.0775259 

0.6250000 

-5001.2847346 

0.8328014 

0.7500000 

-5001.6757710 

0.4417650 

0.8750000 

-5002 . 0145435 

0.1029925 

1.0000000 

-5002 . 1175360 

0.0000000 

i  "Ti"  Vacancy  MIGRATION 

(Min  =  -5113.8368959 

%Distance 

Energy  (eV) 

Difference 

0.0000000 

-5113 . 6320329 

0.2048630 

0.1250000 

-5113.6103820 

0.2265139 

0.2500000 

-5113.5446613 

0.2922346 

0.3750000 

-5113.5031555 

0.3337404 

0.5000000 

-5113.4976717 

0.3392242 

0.6250000 

-5113.5384801 

0.2984158 

0.7500000 

-5113.6303356 

0.2065603 

0.8750000 

-5113.7620650 

0.0748309 

1.0000000 

-5113.8368959 

0.0000000 
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Table  10. 

%Distance 

0.0000000 

0.1250000 

0.2500000 

0.3750000 

0.5000000 

0.6250000 

0.7500000 

0.8750000 

1.0000000 

1.1250000 

1.2500000 

1.3750000 

1.5000000 

1.6250000 

1.7500000 

1.8750000 

2.0000000 

2.1250000 

2.2500000 

2 . 3750000 

2 . 5000000 

2.6250000 

2.7500000 

2.8750000 
3 . 0000000 

3.1250000 

3.2500000 

3 . 3750000 

3.5000000 

3.6250000 

3.7500000 

3 . 8750000 
4.0000000 

4.1250000 

4.2500000 

4.3750000 

4.5000000 

1 . 6250000 

4 .7500000 


FeTi  "Fe"  Vacancy  MIGRATION (Min=-5002 . 5538190) 

Direct  InDirect 

Datal  (eV)  Difference!  Data2  (eV)  Difference2 


-5002.5538190  0.0000000 


-5001.9887234  0.5650956 


-5001.2153024  1.3385166 


-5000.4420590  2.1117600 


-5000.1981781  2.3556409 


-5000.4379935  2.1158255 


-5001.2084583  1.3453607 


-5002.5538190 

0.0000000 

-5002.4249982 

0.1288208 

-5002.0935009 

0.4603181 

-5001.7810563 

0.7727627 

-5001.5660417 

0.9877773 

-5001.4077289 

1.1460901 

-5001.3482255 

1.2055935 

-5001.3337781 

1.2200409 

-5001.1862555 

1.3675635 

-5001.1988910 

1.3549280 

-5000.9026331 

1.6511859 

-5000.8156227 

1.7381963 

-5000.9899902 

1.5638288 

-5001.1942688 

1.3595502 

-5001.5566586 

0.9971604 

-5001.8630311 

0.6907879 

-5001.9657553 

0.5880637 

-5001.8229103 

0.7309087 

-5001.5385912 

1.0152278 

-5001.3129838 

1.2408352 

-5001.2175267 

1.3362923 

-5001.2152742 

1.3385448 

-5001.3265066 

1.2273124 

-5001.4552516 

1.0985674 

-5001.4637261 

1.0900929 

-5001.4271793 

1.1266397 

-5001.2948665 

1.2589525 

-5001.2102149 

1.3436041 

-5001.2176036 

1.3362154 

-5001.2994210 

1.2543980 

-5001.4815633 

1.0722557 

-5001.7413968 

0.8124222 

-5001.9657552 

0.5880638 

-5001.8555345 

0.6982845 

-5001.5487642 

1.0050548 

-5001.1911098 

1.3627092 

-5000.9899979 

1.5638211 

-5000.8166164 

1.7372026 

-5000.8983959 

1.6554231 
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4.8750000 

-5001.1894443 

1.3643747 

5.0000000 

-5001.1864245 

1.3673945 

5.1250000 

-5001.3338043 

1.220014? 

5.2500000 

-5001.9840768 

0.5697422 

-5001.3482314 

1.2055876 

5.3750000 

-5001.4077468 

1.1460722 

5.5000000 

-5001.5660416 

0.9877774 

5.6250000 

-5001.7810170 

0.7728020 

5.7500000 

-5002.0931025 

0.4607165 

5.8750000 

-5002.4248725 

0.1289465 

6.0000000 

-5002.5538138 

0.0000052 

-5002 . 5538185 

0.0000005 
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Table  11.  NiTi  "Ni"  Vacancy  MIGRATION (Min=-5114 . 5049341) 


%Distance 

Direct 

Datal  (eV)  Di 

f ferencel 

InDirect 

Data2  (eV)  Difference2 

0.0000000 

-5114.5049341 

0.0000000 

-5114.5049341 

0.0000000 

0.1250000 

-5114.3898139 

0.1151202 

0.2500000 

-5114.1423992 

0.3625349 

0.3750000 

-5113.9338859 

0.5710482 

0.5000000 

-5113.8152596 

0.6896745 

0.6250000 

-5113.7637272 

0.7412069 

0.7500000 

-5114.1389165 

0.3660176 

-5113.7643917 

0.7405424 

0.8750000 

-5113.7823340 

0.7226001 

1.0000000 

-5113.7839404 

0.7209937 

1.1250000 

-5113.7337982 

0.7711359 

1.2500000 

-5113.6136335 

0.8913006 

1.3750000 

-5113.5094827 

0.9954514 

1.5000000 

-5113 . 5338844 

0.9710497 

-5113.4481337 

1.0568004 

1.6250000 

-5113.4476692 

1.0572649 

1.7500000 

-5113.5085847 

0.9963494 

1.8750000 

-5113.6265550 

0.8783791 

2.0000GCC 

-5113.7095138 

0.7954203 

2.1250000 

-5113.5930761 

0.9118580 

2.2500000 

-5113.0467984 

1.4581357 

-5113 . 3886236 

1.1163105 

2.3750000 

-5113.2578421 

1.2470920 

0.0000000 

-5113.2331188 

1.2718153 

2 . 6250000 

-5113.3017706 

1.2031635 

2.7500000 

-5113.4426337 

1.0623004 

2.8750000 

-5113.5775667 

0.9273674 

3 . 0000000 

-5112.8455693 

1.6593648 

-5113.5894853 

0.9154488 

3 . 1250000 

-5113 . 5508028 

0.9541313 

3.2500000 

-5113.4123014 

1.0926327 

3.3750000 

-5113.2914537 

1.2134804 

3.5000000 

-5113.2331181 

1.2718160 

3 . 6250000 

-5113.2515763 

1.2533578 

3.7500000 

-5113 . 0467884 

1.4581457 

-5113.3547116 

1.1502225 

3.8750000 

-5113.5363502 

0.9685839 

4 . 0000000 

-5113.7095137 

0.7954204 

4.1250000 

-5113.6438446 

0.8610895 

4.2500000 

-5113.5184936 

0.9864405 

4.3750000 

-5113.4489913 

1.0559428 

4 . 5000000 

-5113.5338719 

0.9710622 

-5113.4481258 

1 . 0568083 

4.6250000 

-5113.5134446 

0 . 9914895 

4.7500000 

-5113.6249491 

0.8799850 
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4.8750000 

-5113.7477682 

0.7571659 

5.0000000 

-5113.7839404 

0.7209937 

5.1250000 

-5113.7823341 

0.7226000 

5.2500000 

-5114.1389188 

0.3660153 

-5113.7643917 

0.7405424 

5.3750000 

-5113.7637272 

0.7412069 

5.5000000 

-5113.8152596 

0.6896745 

5.6250000 

-5113.9388490 

0.5660851 

5.7500000 

-5114.1423971 

0.3625370 

5.8750000 

-5114.3898138 

0.1151203 

6.0000000 

-5114.5049341 

0.0000000 

-5114 . 5049341 

0.0000000 
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Table  12.  Stable  Defect  Configuration  Energies  Results 

(E  =  U  -  Uo) 

Nickel  Titanium  (NiTi) 


System 

Energy (eV) 

Difference (eV) 

Perfect  Crystal  (Uo) 

-5120.8607297 

0.0000000 

Ni  Antisite  Defect 

-5120.0983898 

0.7623399 

Ti  Antisite  Defect 

-5120.6076298 

0.2530999 

NiTi  Antisite  Formation 

0.5092400 

NiTi  Antisite  Pair 

1.0154398 

Ni  Vacancy  Defect 

-5114.5049341 

6.3557956 

Ti  Vacancy  Defect 

-5113 .8368959 

7.0238338 

NiTi  Vacancy  Formation 

1.4800000 

Ni  Interstitial  Defect 

-5123.2559795 

-2.3952498 

Ti  Interstitial  Defect 

-5123.4626653 

-2 . 6019356 

Ni  Frenkel  Pair  in  NiTi 

3.9605458 

Ti  Frenkel  Pair  in  NiTi 

4.4218982 

Iron  Titanium  (FeTi) 


System 

Energy (eV) 

Difference (eV) 

(U) 

(L) 

Perfect  Crystal  (Uo) 

-5008 .4530840 

0.0000000 

Fe  Antisite  Defect 

-5007.9620535 

0.4910305 

Ti  Antisite  Defect 

-5008 .2178813 

0.2352027 

FeTi  Antisite  Formation 

0.2558278 

FeTi  Antisite  Pair 

0.7262332 

Fe  Vacancy  Defect 

-5002 . 5538190 

5.8992650 

Ti  Vacancy  Defect 

-5002 . 1175360 

6.3355480 

FeTi  Vacancy  Formation 

1.0700000 

Fe  Interstitial  Defect 

-5010.0113699 

-1.5582859 

Ti  Interstitial  Defect 

-5010.3980165 

-1.9449325 

Fe  Frenkel  Pair  in  FeTi 

4 .3409791 

Ti  Frenkel  Pair  in  FeTi 

4.3906155 
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APPENDIX  C: 


Vacancy  Formation  Energy  Code 


The  following  program  was  written  using  TK  Solver  Plus 
software  version  1.1.  Eight  unknowns  NTlt  n J',  n „l,  nral , 

nTal ,  and  | iT‘  are  solved  by  the  eight  equations  (8a,  8b), 

(9a, 9b),  (10a, 10b),  (11),  (12),  and  their  complements  as 
presented  in  section  V.  A  list  of  guessed  values  are  input 
into  the  program  for  each  of  the  unknowns  for  each  tempera¬ 
ture  of  interest.  The  only  required  input  values  to  the 
program  are  the  five  total  system  energies  for  the  global 
minimum  configurations  and  the  number  of  lattice  sites  for 
each  type  of  atom  in  the  compound.  These  energies  include 
the  perfect  crystal  U0  for  each  alloy,  the  antisite  defect 
total  system  energies  for  NiTi  b'al,UTal  or  for  FeTi  UFae  ,UTal 
(Table  1) ,  plus  the  vacancy  defect  total  system  energies 
for  NiTi  Uf'.Ul'  or  for  FeTi  U[e  ,UTvl  (Table  3). 

The  program  is  written  using  nomenclature  of  the  NiTi 
lattice.  This  same  approach  was  used  in  the  presentation 
of  the  equations  in  section  V.  The  program  is  however 
generic  and  can  be  used  for  any  type  of  compound. 
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S  S  £  S  R  K 


======  VARIABLE  SHEET  ===============  For  Academic  Use  Only 

St  Input  Name  Output  Unit  Comment======================== 

•k’k'klc'k’k'k’k'k'kick'k'k'k'k'k'kit'k'k'k'k'kit’k’k'k'k'k'k'k 

TK  Solver  Program  Ver5.4  12Dec90 
******************************** 

Calculates  the  effective  vacancy 
formation  energy,  using  Foiles 
&  Daw  method (J. Mat .Res . 2 , 3 , 1987) 
File:  ENERGY. TK  by  Tom  button 

******************************** 


L 

L 


LG 


m 

▲ 

1200  K 

Temperature  (Kelvin) 

.83333 

TInv 

1/K 

Temperature  Inverse  =  1000/T 

OUTPUT: 

Cv 

3.0E-7 

Concentration  of  Vacancy  Defect 
GUESS  INPUT  VALUES: 

2.9E-7 

n  Ni  v 

Atom  fraction  Ni  Vacancy  Defect 

5 . 5E-9 

n  Ti  v 

Atom  fraction  Ti  Vacancy  Defect 

. 00366 

n  Ni  a 

Atom  fraction  Ni  Antisite  Defect 

.00366 

n  Ti  a 

Atom  fraction  Ti  Antisite  Defect 

-4.874 

muNi 

eV 

Chemical  Potential  of . Ni 

-5.128 

muTi 

eV 

Chemical  Potential  of . Ti 

511.99 

N_Ni 

#  Lattice  Sites (@given  T) Ni 

511.99 

NJTi 

#  Lattice  Sites (@given  T) Ti 
INPUT : 

Chart  of  Nuclides (k=8. 61735E-5) 

8.6E-5 

k 

eV/K 

Boltzmann  Constant  (eV/Kelvin) 

1024 

Ns 

Total  #  Lattice  Sites (@  T=0) 

512 

nsNi 

#  Lattice  Sites (@  T=0)  of  Ni 

512 

nsTi 

#  Lattice  Sites (@  T=0)  of  Ti 

fsNi 

.5 

Fraction  of  Ni  Lattice  Sites 

fsTi 

.5 

Fraction  of  Ti  Lattice  Sites 

i 

2 

Inverse  Fraction  =  1/fsNi 

j 

2 

Inverse  Fraction  =  1/fsTi 
INPUT  ENERGIES: 

-5120.9 

Uo 

eV 

NiTi  System  Energy  Perfect  Cryst 

-5120.1 

U  Ni  a 

eV 

NiTi  System  Energy  Ni  Antisite 

-5120.6 

U  Ti  a 

eV 

NiTi  System  Energy  Ti  Antisite 

-5114 . 5 

U_N i_v 

eV 

NiTi  System  Energy  Ni  Vacancy 

-5113.8 

U_Ti_v 

eV 

NiTi  System  Energy  Ti  Vacancy 
OUTPUTS : 

E  Ni  a 

.7623  eV 

NiTi ...  Energy  of  Ni  Antisite  Def 
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E_Ti_a  .2530  eV  NiTi ...  Energy  of  Ti  Antisite  Def 

E_Ni_v  5.355  eV  NiTi ...  Energy  of  Ni  Vacancy  Def 

E_Ti_v  7.02  3  eV  NiTi ...  Energy  of  Ti  Vacancy  Def 

U  -5117.  eV  NiTi  Total  System  Energy  Ns  Site 

S  44.32  eV/K  Configuration  Entropy  NiTi  Syste 

Ideal  entropy  Function 
si  .0432  eV/K  Term  #1 

s2  9.1E-6  eV/K  Term  #2 

s3  .0433  eV/K  Term  #3 

s4  2.1E-7  eV/K  Term  #4 

OUTPUT : 

Vacancy 

L  Efxslop  Energy  Formation=-k*ln(Cv)/TA-l 

NOTE:  The  following  values 

are  output  initially&later 
used  as  inputs  (see  Rules) . 

L  TInvl  .8333  Dummy  List  Offset  TInv  by  1 

L  Cvl  3.0E-7  Dummy  List  Offset  Cv  by  1 
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************************************* 


-5120 

-5120 

-5120 

-5114 

-5113 


-5008 

-5007 

-5008 

-5002 

-5002 


-4404 

-4403 

-4405 

-4399 

-4398 


-3965 

-3964 

-3966 

-3960 

-3959 


INPUT  VALUES: 

************************************* 


8607 

NiTi 

System 

Energy. 

.Perfect  Crystal 

0984 

NiTi 

System 

Energy. 

.Ni 

Antisite 

Defect 

6076 

NiTi 

System 

Energy. 

.Ti 

Antisite 

Defect 

5045 

NiTi 

System 

Energy. 

.Ni 

Vacancy 

Defect 

8369 

NiTi 

System 

Energy. 

•  Ti 

Vacancy 

Defect 

Note : 

Ns=1024  nsNi=512  nsTi=512 

4531 

FeTi 

System 

Energy. 

.Perfect  Crystal 

9621 

FeTi 

System 

Energy . 

.Ni 

Antisite 

Defect 

2179 

FeTi 

System 

Energy. 

.Ti 

Antisite 

Defect 

5538 

FeTi 

System 

Energy. 

•  Ni 

Vacancy 

Defect 

1175 

FeTi 

System 

Energy. 

.Ti 

Vacancy 

Defect 

Note : 

Ns=1024  nsFe=512  nsTi=512 

4207 

CuTi 

System 

Energy. 

.Perfect  Crystal 

0266 

CuTi 

System 

Energy. 

.Ni 

Antisite 

Defect 

4299 

CuTi 

System 

Energy. 

.Ti 

Antisite 

Defect 

6440 

CuTi 

System 

Energy. 

.Ni 

Vacancy 

De.  ~ct 

4100 

CuTi 

System 

Energy. 

.Ti 

Vacancy 

Defect 

Note:  Ref  Case  used  k=8.615E-5 
Note:  Ns=1024  nsCu=512  nsTi=512 


5986 

CuTi  2 

System 

Energy 

Perfect  Crystal 

3864 

CuTi  2 

System 

Energy 

Ni 

Antisite 

Defect 

3511 

CuTi2 

System 

Energy 

Ti 

Antisite 

Defect 

8716 

CuTi  2 

System 

Energy 

Ni 

Vacancy 

Defect 

7368 

CuTi  2 

System 

Energy 

Ti 

Vacancy 

Defect 

Note : 

Ref  Case 

used  k=8 

.  615E-5 

Note:  Ns=882  nsCu=294  nsTi=588 
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========  RULE  SHEET  =================  For  Academic  Use  Only 

^ ********************************************************** 
"Equations  for  NiTi  Effective  Vacancy  Formation  Energy 

"A********************************************************* 

********************* *TOTAL  VACANCY  CONCENTRATION:********** 
Cv  =  n_Ni_v  +  n_Ti_v 

'•*********************** INVERSE  TEMPERATURE :*************** 
TInv  =  1000  /  T 

•«*****INPUT  (Total  System  Energy )- (Perfect  Crystal  Energy) 

E_Ti_a  =  U_Ti_a  -  Uo 

E_N i_a  =  U_N i_a  -  Uo 

E_Ti_v  =  U_Ti_v  -  Uo 

E_Ni_v  =  U_Ni_v  -  Uo 

»*****jnPUT  of  Number  of  Lattice  Sites  of  Each  Element:**** 
fsNi  =  nsNi  /  (nsNi  +  nsTi) 
fsTi  =  nsTi  /  (nsNi  +  nsTi) 

***************************  *  *EQUATION . .  (1) ***************** 
n_Ni_v=  fsNi* (exp (- (E_Ni_v+muNi)/ (k*T) ) )  / 

( 1+exp (- ( E_N i_v+muN i ) / (k*T) ) ) 

"********************* *******EQUATION. . (2)***************** 
n_Ti_v=  fsTi* (exp (- (E_Ti_v+muTi)/ (k*T) ) )  / 

(1+exp (- (E_Ti_v+muTi)/ (k*T) ) ) 

11  ***************************  *  EQUATION .  .  (3)***************** 
n _Ni_a=fsTi* (exp (- (E_Ni_a+muTi-muNi) / (k*T) ) )  / 

( 1+exp (- (E_Ni_a+muTi-muNi)  /  (k*T) ) ) 

*>  ****************************  EQUATION .  .  (4)  ***************** 
n  Ti_a=fsNi* (exp(- (E_Ti_a+muNi-muTi) / (k*T) ) )  / 

(1+exp (- (E_Ti_a+muNi-muTi)  /  (k*T) ) ) 

******** *********************EQUATI0N. . (5)***************** 

N_Ni  =  Ns  *  (fsNi  -  n_Ni_v  +  n_Ni_a  -  n_Ti_a) 

**************************  *  *EQUATION  .  .  (6)  ***************** 
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N_Ti  =  Ns  *  ( fsTi  -  n_Ti_v  +  n_Ti_a  -  n_Ni_a) 

•» ************************* ***EQUATION.  .  (7)***************** 
fsNi  =  N_Ni  /  (N_Ni  +  N_Ti) 

****************************  *  *  EQUATION .  .  (8)  ***************** 
U  -  (k*T*S)  =  (muNi*N_Ni)  +  (muTi*N_Ti) 


"WHERE : 

U  =U  o+N  s  * ( n_N i_v  *  E_N i_v  +  n_Ti_v*E_Ti_v  +  n_Ni_a*E_Ni_a  + 

n_Ti_a*E_Ti_a) 

S  =Ns  *  ( (fsNi) * (sl+s2)  +  (fsTi) * (s3+s4) ) 

i=l/fsNi  "Temporary  Term 
j=l/fsTi  "Temporary  Term 
If  n_Ti_a>0  then  sl=  -( (i*n_Ti_a) *ln(i*n_Ti_a)+ 

(l-i*n_Ti_a)  *  In (l-i*n_Ti_a) ) 

*If  n_Ti_a<=0  then  sl=  0 

If  n_Ni_v>0  then  s2=  - ( ( i*n_Ni_v) *ln ( i*n_Ni_v) + 

(l-i*n_Ni_v)  *  In ( l-i*n_Ni_v) ) 

*If  n_Ni_v<=0  then  s2=  0 

If  n_Ni_a>0  then  s3=  - ( ( j *n_Ni_a) *ln ( j *n_Ni_a) + 

(l-j*n_Ni_a)  *  ln(l-j*n_Ni_a) ) 

*If  n_Ni_a<=0  then  s3=  0 

If  n_Ti_v>0  then  s4=  - ( ( j*n_Ti_v) *ln( j *n_Ti_v) + 

( 1- j  *n_Ti_v) *ln ( 1- j  *n_Ti_v) ) 

*If  n  Ti  v<=0  then  s4=  0 


"Note:  The  If  statements  were  added  to  help  cure  natural 

log  errors,  which  occur  when  the  atom  fraction  values  go 
negative  during  iteration.  The  negative  values  occur  when 
errors  become  larger  than  the  tolerances  of  the  iterations. 
Final  values  must  be  positive.  This  is  not  a  solution  to 
the  problem  of  log  errors,  but  a  symptom  of  low  tolerances. 
Solutions  include  changing  the  input  guess  values,  changing 
the  tolerance  (default=lE-6) ,  and  changing  the  form  of  the 
equations.  In  general,  changing  the  equation  doesn't 
effect  results.  The  first  solution  is  not  always  correct; 
change  tolerances  and  test  the  solution.  Problems  with 
"too  many  guesses",  negative  atom  fraction  values,  incon¬ 
sistency,  and  more  than  one  solution  are  all  symptoms  of 
low  tolerances  compared  to  errors.  For  understanding  the 
problems  read  the  TK  Solver  help  sections  on  the  iterative 
solver . 
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The  problems  are  mostly  with  the  antisite  defect  atom 
fractions  and  they  become  more  pronounced  with  lower  tem¬ 
peratures.  Thus,  to  solve  these  equations  it  is  best  to 
start  with  a  higher  temperature  input  (1200  Kelvin)  and 
find  a  solution.  Use  the  guess  values  given  in  the  vari¬ 
able  sheet  to  start  with.  Then  use  the  results  as  input  to 
the  next  lower  temperature  of  interest.  This  ensures 
reasonable  input  values.  As  solution  errors  appear  lower 
the  tolerance.  At  some  point  the  errors  are  so  small  com¬ 
pared  to  the  tolerances,  that  no  solution  is  possible  with 
TK.  At  these  low  temperatures,  the  values  are  generally 
not  realistic  and  the  process  can  be  stopped. 

,,***************************^t****************************** 

"  *For  List  Solutions  VACANCY  FORMATION  ENERGY (Ef)=  -k*slope 

•*  a********************************************************* 

"Note:  First  calculate  Cv  &  TInv  valves.  Then  offset  these 
values  by  1  space  value  at  top  of  each  list.  Next  cancel 
the  TInv=  &  Cv=  rules  &  uncancel  the  Efxslope  rule  to  do 
the  calculation.  Be  sure  to  change  Cv,  Cvl,  TInvl,  &  TInv 
to  list  Inputs. 


C  TInv  =  TInvl 
Cv  =  Cvl 

C  Efxslope=  -k  *  (In (Cvl) -ln(Cv) )  /  ( (Tlnvl-TInv) /1000) 


=====  UNIT  SHEET  ============  For  Academic  Use  Only 

To=======  Multiply  By===Add  Of f set===Comment===== 

C  -273.15 

F  1.8  32 

F  1.8  -459.67 

eV*KA-l 
eV 


From===== 

K 
C 
K 

eV/K 
eV/K*K 

=============  GLOBAL  sheet  = 

Display  Intermediate  Values: 
Stop  on  List  Error: 

Use  Automatic  Iteration: 
Comparison  Tolerance: 

Typical  Value: 

Maximum  Iteration  Count: 


========  For  Academic  Use  Only 

Yes 

No 

Yes 

IE-10 

1 

25 
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T(K) 

10000 

5000 

3333 

2500 

2000 

1667 

1429 

1250 

1200 

1111 

1000 

909.1 

833.3 
750 

714.3 
666.7 

571.4 

545.5 
500 
444.4 

428.6 
400 
363  .  S 
352.9 
333.3 

307.7 
300 

285.7 

Key: 

T 

KA-1 

CvNiTi 

n_Ni_' 

n_Ni_j 

muNi 

NNi 


Table  13.  Vacancy  Concentration  of  NiTi 


KA- 

1  n_N i _v 

n_Ti_v 

n_Ni_a 

r._Ti_a  rouNi  muTi 

NNi 

.  1 

. 12827163 

.09360703 

.18717 

.16984  -5.43  -5.75 

398 

.2 

.02036074 

.00826969 

. 12072 

.11468  -4.99  -5.26 

497 

.  3 

. 00333461 

.00080280 

. 07356 

.07229  -4.91  -5.17 

510 

.4 

. 00056144 

8 . 2722E-5 

.04339 

.04315  -4.89  -5.14 

512 

.5 

9 . 6792E-5 

8.8014E-6 

.02499 

.02494  -4.88  -5.13 

512 

.6 

1. 6952E-5 

9 . 5332E-7 

. 01417 

.01416  -4.87  -5.13 

512 

.7 

2 . 9973E-6 

1 . 0430E-7 

.00796 

.00796  -4.87  -5.13 

512 

.8 

5 . 3287E-7 

1 . 1477E-8 

. 00445 

.00445  -4.87  -5.12 

512 

.83 

2 . 9987E-7 

5 . 5038E-9 

.00366 

.00366  -4.87  -5.12 

512 

.9 

9 . 5034E-8 

1 . 2  668E-9 

.00248 

.00248  -4.87  -5.12 

512 

1 

1.6978E-8 

1 . 401E-10 

. 00138 

.00138  -4.87  -5.12 

512 

1.1 

3 . 0362E-9 

1 . 550E-11 

.00076 

.00076  -4.87  -5.12 

512 

1.2 

5 . 433E-10 

1 . 717E-12 

.00042 

.00042  -4.87  -5.12 

512 

1.3 

5 . 480E-11 

9. 136E-14 

.00019 

.00019  -4.87  -5.12 

512 

1.4 

1 . 741E-11 

2 . 107E-14 

. 00013 

.00013  -4.87  -5.12 

512 

1.5 

3.116E-12 

2 . 335E-15 

7.3E-5 

7.3E-5  -4.87  -5.12 

512 

1.7 

4 . 227E-14 

9 . 546E-18 

1.7E-5 

1.7E-5  -4.87  -5.12 

512 

1.8 

1. 008E-14 

1.526E-18 

.00001 

.00001  -4.87  -5.12 

512 

2 

5.734E-16 

4 . 033E-20 

3.8E-6 

3.8E-6  -4.87  -5.12 

512 

2.2 

7 . 7  65E-18 

1 . 078E-20 

8.7E-7 

8.7E-7  -4.87  -5.12 

512 

2.3 

1.751E-18 

2 . 981E-20 

5.4E-7 

5.4E-7  -4.87  -5.12 

512 

2.5 

6 . 838E-20 

3 . 542E-20 

2.0E-7 

2.0E-7  -4.87  -5.12 

512 

2.7 

1. 166E-21 

2 . 838E-22 

4.6E-8 

4.6E-8  -4.87  -5.12 

512 

2.8 

2.612E-21 

-2.27E-21 

2.8E-8 

2.8E-8  -4.87  -5.12 

512 

3 

-1.81E-19 

2 . 205E-20 

1. IE-8 

1.1E-8  -4.87  -5.12 

512 

3 . 2 

-3.55E-17 

-2 . 95E-19 

2.4E-9 

2.4E-9  -4.87  -5.12 

512 

3.3 

-6.21E-17 

-2.18E-20 

1.5E-9 

1.5E-9  -4.87  -5.12 

512 

3.5 

3 . 371E-16 

7. 154E-19 

6E-10 

6E-10  -4.87  -5.12 

512 

Temperature  in  Kelvin 
=  Inverse  Temperature  *  1000  in  1  /Kelvin 
=  Total  Vacancy  Concentration 

/  or  n_Ti_v  =  Ni  or  Ti  Vacancy  Concentration  Fraction 
i  or  n_Ti_a  =  Ni  or  Ti  Antisite  Defect  Concentration  Fraction 
or  muTi  =  Ni  or  Ti  Chemical  Potential 
=  Ni  or  Ti  Number  of  lattice  sites  at  a  given  Temperature  (T) 
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Table  14.  Vacancy  Concentration  of  FeTi 


T  (K) 

KA- 

1  n_Fe_v 

n_Ti_v 

n_Fe_a 

n_Ti_a  muFe  muTi 

NFe 

10000 

.  1 

.21754061 

. 18324654 

.20673 

.18958  -5.67  -5.86 

307 

5000 

.2 

.  05608482 

.03077715 

.15688 

.14423  -5.00  -5.16 

468 

3333 

.3 

. 01498607 

.00537199 

.11255 

.10775  -4.90  -5.03 

502 

2500 

.4 

.00405130 

.00098554 

.07895 

.07742  -4.86  -4.99 

509 

2000 

.5 

. 00111089 

.00018685 

. 05445 

.05398  -4.84  -4.97 

511 

1667 

.6 

.00030883 

.00003617 

.03702 

.03688  -4.83  -4.96 

512 

1429 

.7 

8 . 6802E-5 

7 . 0972E-6 

.02490 

.02486  -4.83  -4.96 

512 

1250 

.8 

2 . 4594E-5 

1 . 4051E-6 

.01661 

.01660  -4.83  -4.95 

512 

1200 

.83 

1 . 6175E-5 

8 . 2012E-7 

. 01450 

.01449  -4.83  -4.95 

512 

1111 

.9 

7 . 0072E-6 

2 . 7985E-7 

.01102 

.01102  -4.82  -4.95 

512 

1000 

1 

2 . 0041E-6 

5 . 5955E-8 

. 00729 

.00729  -4.82  -4.95 

512 

909.1 

1.1 

5.7465E-7 

1. 1218E-8 

.00481 

.00481  -4.82  -4.95 

512 

833.3 

1.2 

1 . 6506E-7 

2 . 2527E-9 

. 00316 

.00316  -4.82  -4.95 

512 

750 

1.3 

3 . 1333E-8 

2 . 654E-10 

.00181 

.00181  -4.82  -4.95 

512 

714.3 

1.4 

1.3658E-8 

9. 112E-11 

.00137 

.00137  -4.82  -4.95 

512 

666.7 

1.5 

3.9322E-9 

1.834E-11 

.00090 

.00090  -4.82  -4.95 

512 

571.4 

1.7 

1.751E-10 

3. 339E-13 

.00031 

.00031  -4.82  -4.95 

512 

545.5 

1.8 

6.207E-11 

8 . 784E-14 

.00022 

.00022  -4.82  -4.95 

512 

500 

2 

7 . 803E-12 

6 . 082E-15 

.00011 

.00011  -4.82  -4.95 

512 

444.4 

2.2 

3 . 478E-13 

1 . 108E-16 

3 . 8E-5 

3.8E-5  -4.82  -4.95 

512 

428.6 

2.3 

1 . 233E-13 

2 . 916E-17 

2 . 7E-5 

2.7E-5  -4.82  -4.95 

512 

400 

2.5 

1 . 550E-14 

2 . 019E-18 

1.3E-5 

1.3E-5  -4.82  -4.95 

512 

363.6 

2.7 

6 . 912E-16 

3 . 681E-20 

4.6E-6 

4.6E-6  -4.82  -4.95 

512 

352.9 

2.8 

2 . 451E-16 

9 . 662E-2 1 

3. 3E-6 

3.3E-6  -4.82  -4.95 

512 

333.3 

3 

3 . 081E-17 

4 . 874E-22 

1. 6E-6 

1.6E-6  -4.82  -4.95 

512 

307.7 

3.2 

1 . 4 11E-18 

3 . 762E-20 

5. 6E-7 

5.6E-7  -4.82  -4.95 

512 

300 

3.3 

4 . 578E-19 

1 . 788E-21 

4.0E-7 

4. OE-7  -4.82  -4.95 

512 

285.7 

3.5 

1. 626E-19 

-9.35E-20 

2 . OE-7 

2. OE-7  -4.82  -4.95 

512 

Key: 

T  =  Temperature  in  Kelvin 

KA-1  =  Inverse  Temperature  *  1000  in  1  /Kelvin 

CvFeTi  =  Total  Vacancy  Concentration 

n_Fe_v  or  n_Ti_v  =  Fe  or  Ti  Vacancy  Concentration  Fraction 
n_Fe  a  or  n__Ti_a  =  Fe  or  Ti  Antisite  Defect  Concentration  Fraction 
muFe  or  muTi  =  Fe  or  Ti  Chemical  Potential 
NFe  =  Fe  or  Ti  Number  of  lattice  sites  at  a  given  Temperature  (T) 
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Table  15.  Vacancy  Concentration  of  CuTi 
T(K)  KA-1  n_Cu_v  n_Ti_v  n_Cu_a  n_Ti_a  muCu  muTi  NCu 


10000 

.  1 

5000 

.2 

3333 

.3 

2500 

.4 

2000 

.5 

1667 

.  6 

1429 

.7 

1250 

.8 

1200 

.83 

1111 

.9 

1000 

1 

909.1 

1.1 

833 . 3 

1.2 

750 

1.3 

714.3 

1.4 

666.7 

1.5 

571.4 

1.7 

545.5 

1.8 

500 

2 

444.4 

2.2 

428.6 

2.3 

400 

2.5 

363.6 

2.7 

352.9 

2.8 

333 . 3 

3 

307.7 

3.2 

300 

3.3 

285.7 

3.5 

.25516498 
. 06687083 
.01776014 
.00475944 
. 00128679 
. 00035093 
9 . 6472E-5 
2 . 6709E-5 
1 . 7432E-5 
7 . 4397E-6 
2 . 0831E-6 
5 . 8583E-7 
1 . 6535E-7 
3 . 0755E-8 
1 . 3285E-8 
3.7778E-9 
1.640E-10 
5.772E-11 
7 . 161E-12 
3 . 139E-13 
1. 107E-13 
1 . 379E-14 
6.062E-16 
2.140E-16 
2.667E-17 
1. 174E-18 
4. 146E-19 
5 . 125E-20 


.25142088 
.06307794 
. 01598982 
.00410771 
.00106767 
. 00028025 
7 . 4187E-5 
1. 9781E-5 
..00001275 
5. 3068E-6 
1.4311E-6 
3 . 8765E-7 
1. 0538E-7 
1 . 8644E-8 
7.8543E-9 
2 . 1512E-9 
8 . 501E-11 
2.900E-11 
3. 380E-12 
1 . 349E-13 
4 . 611E-14 
5 . 393E-15 
2 . 159E-16 
7 . 387E-17 
8 . 646E-18 
3 . 466E-19 
1 . 186E-19 
1 . 353E-20 


.22314 
. 19603 
.16969 
.14538 
.12337 
.  10378 
.08659 
.07175 
.06729 
.  05908 
.04840 
.03948 
. 03208 
.  02422 
. 02101 
.01695 
.00984 
.00820 
. 00568 
.00326 
.  00271 
.00187 
.00107 
.00089 
.00061 
.00035 
.  00029 
.  00020 


.22126 
.19413 
.  16881 
.14505 
.12326 
.10374 
. 08658 
.07174 
.  06729 
. 05908 
.04840 
.  03948 
. 03208 
. 02422 
.02101 
.01695 
.00984 
.00820 
.  00568 
.00326 
.  00271 
.00187 
.  00107 
.00089 
. 00061 
.  00035 
.  00029 
.  00020 


-4.81 

-3.97 

-3.82 

-3.77 

-3.75 

-3.73 

-3.72 

-3.71 

-3.71 

-3.71 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 

-3.70 


-6.02 

253 

-5.17 

445 

-5.03 

495 

-4.97 

507 

-4.95 

511 

-4.93 

512 

-4.92 

512 

-4.91 

512 

l 

.t* 

• 

VO 

i-1 

512 

-4.91 

512 

-4.91 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

-4 . 90 

512 

-4.90 

512 

-4.90 

512 

-4.90 

512 

Key: 

T  =  Temperature  in  Kelvin 

KA-1  =  Inverse  Temperature  *  1000  in  1/Kelvin 

CvCuTi  =  Total  Vacancy  Concentration 

n_Cu__v  or  n_Ti_v  =  Cu  or  Ti  Vacancy  Concentration  Fraction 
n_Cu_a  or  n_Ti_a  =  Cu  or  Ti  Antisite  Defect  Concentration  Fraction 
muCu  or  muTi  =  Cu  or  Ti  Chemical  Potential 
NCu  =  Cu  or  Ti  Number  of  lattice  sites  at  a  given  Temperature  (T) 
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Table  16.  Vacancy  Concentration  CuTi2 


T  (K) 

KA-1 

n_Cu_v 

n_Ti_v  n 

_Cu_a 

n_Ti_a  muCu 

muTi2 

NCu 

10000 

.  1 

.21562640 

.32798014 

.21719 

.18277  -5.24 

-5.83 

134 

5000 

.2 

.07974555 

.08346786 

.18367 

.15833  -4.22 

-5.02 

246 

3333 

.3 

.02803006 

.02225664 

.14880 

.13753  -4.04 

-4.89 

279 

2500 

.4 

. 00954481 

.00616024 

.12062 

.11631  -3.96 

-4.85 

289 

2000 

.5 

, 00322218 

.00174452 

.09762 

.09605  -3.92 

-4.83 

293 

1667 

.6 

.00109218 

.00050101 

.07850 

.07794  -3.90 

-4.82 

294 

1429 

.7 

.00037326 

.00014527 

.06263 

.06243  -3.89 

-4.82 

294 

1250 

.8 

. 00012866 

4 . 2427E-5 

.04958 

.04951  -3.88 

-4.82 

294 

1200 

.83 

9 . 0371E-5 

2 . 8185E-5 

.04579 

.04574  -3.87 

-4.82 

294 

1111 

.9 

4 . 4689E-5 

.00001246 

.03899 

.03896  -3.87 

-4 .81 

294 

1000 

1 

1 . 5621E-5 

3 . 6756E-6 

.03048 

.03048  -3.86 

-4.81 

294 

909.1 

1.1 

5 . 4880E-6 

1 . 0881E-6 

.02373 

.02373  -3.86 

-4.81 

294 

833 . 3 

1.2 

1 . 9360E-6 

3 . 2300E-7 

.01840 

.01840  -3.86 

-4.81 

294 

750 

1.3 

4.8489E-7 

6 . 4175E-8 

.01306 

.01306  -3.85 

-4.81 

294 

714.3 

1.4 

2.4306E-7 

2 . 8637E-8 

.01098 

.01098  -3.85 

-4.81 

294 

666.7 

1.5 

8.6398E-8 

8 . 5456E-9 

.00846 

.00846  -3.85 

-4.81 

294 

571.4 

1.7 

6 . 5463E-9 

4. 174E-10 

.00438 

.00438  -3.85 

-4.81 

294 

545.5 

1.8 

2.7736E-9 

1 . 527E-10 

.00352 

.00352  -3.85 

-4.81 

294 

500 

2 

4 . 985E-10 

2.046E-11 

.00226 

.00226  -3.85 

-4 .81 

294 

444.4 

2.2 

3 . 806E-11 

1 . 004E-12 

.00116 

.00116  -3.85 

-4 .81 

294 

428.6 

2.3 

1. 615E-11 

3 . 679E-13 

.00093 

.00093  -3.85 

-4.81 

294 

400 

2.5 

2 . 910E-12 

4 . 936E-14 

. 00060 

.00060  -3.84 

-4.82 

294 

363.6 

2.7 

2.226E-13 

2 . 427E-15 

.00031 

.00031  -3.84 

-4 . 82 

294 

352.9 

2.8 

9 . 451E-14 

8 . 892E-16 

.00025 

.00025  -3.84 

-4.82 

294 

333.3 

3 

1.704E-14 

1. 194E-16 

. 00016 

.00016  -3.84 

-4.82 

294 

307.7 

3.2 

1. 304E-15 

5 . 872E-18 

8. IE-5 

8.1E-5  -3.84 

-4.82 

294 

300 

3 . 3 

5. 537E-16 

2 . 151E-18 

6. 5E-5 

6.5E-5  -3.84 

-4.82 

294 

285.7 

3.5 

9 . 983E-17 

2 . 888E-19 

4 . 2E-5 

4.2E-5  -3.84 

-4.82 

294 

Key: 

T  =  Temperature  in  Kelvin 

KA-1  =  Inverse  Temperature  *  1000  in  1 /Kelvin 

CvCuTi2  =  Total  Vacancy  Concentration 

n_Cu_v  or  n_li_v  =  Cu  or  Ti2  Vacancy  Concentration  Fraction 
n_Cu_a  or  n_Ti_a  =  Cu  or  Ti2  Antisite  Defect  Concentration  Fraction 
muCu  or  muTi2  =  Cu  or  Ti2  Chemical  Potential 

NCu  =  Cu  or  Cu*2  =  Ti2  Number  of  lattice  sites  at  a  given  Temperature  (T) 
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Table  17.  Vacancy  Concentration  of  FeTi,NiTi,CuTi, &CuTi2 


Temp(K)  KA-1  CvCuTi2  CvCuTi  CvFeTi  CvNiTi 


10000 

.  1 

5000 

.2 

3333.3 

.3 

2500 

.4 

2000 

.  5 

1666.7 

.6 

1428.6 

.7 

1250 

.8 

1200 

.833 

1111.1 

.9 

1000 

1 

909.1 

1.1 

833.3 

1.2 

750 

1.33 

714.3 

1.4 

666.7 

1.5 

571.4 

1.75 

545.5 

1.83 

500 

2 

444.4 

2 .25 

428.6 

2.33 

400 

2.5 

363.6 

2.75 

352.9 

2.83 

333.3 

3 

307.7 

3.25 

300 

3.33 

285.7 

3.5 

. 5436065413 
.1632134099 
.0502867013 
.0157050543 
.0049667039 
.0015931990 
. 0005185351 
.0001710892 
.0001185558 
5 . 7 14868E-5 
1. 929607E-5 
6 . 576084E-6 
2 . 258967E-6 
5 . 490644E-7 
2 . 716983E— 7 
9 . 494356E-8 
6. 963697E-9 
2 . 926309E— 9 
5. 18962E-10 
3 . 90646E— 11 
1 . 65196E-11 
2 . 95910E-12 
2 . 25032E-13 
9 . 53972E-14 
1 . 7 1555E-14 
1 . 30990E-15 
5 . 55843E-16 
1 . 00114E-16 


.5065858619 
.1299487611 
.0337499548 
.0088671528 
.0023544524 
.0006311847 
.0001706591 
4 . 648931E-5 
3 . 018146E-5 
1. 274647E-5 
3 . 514274E-6 
9 . 734784E-7 
2 . 7  07  338E-7 
4 . 939816E-*5 
2 . 113923E-8 
5 . 929085E-9 
2.48986E-10 
8.67145E-11 
1. 05413E-11 
4 . 48752E-13 
1 . 5683  6E-13 
1 . 91788E-14 
8 . 22 104E-16 
2 . 87  87  4E-16 
3 . 53 17  3E-17 
1 . 52076E-18 
5. 33225E-19 
6 . 477  63E-2  0 


.4007871513 
.0868619735 
.0203580644 
.0050368411 
. 0012977370 
.0003449970 
.0000938990 
2 . 599875E-5 
1. 699486E-5 
7 . 287051E-6 
2 . 060060E-6 
5 . 8  58714E-7 
1. 673087E-7 
•>  .159793E-8 
1 . 374903E— 8 
3 . 950536E-9 
1.75429E-10 
6 . 2 1618E-11 
7 . 80877E-12 
3 . 47911E-13 
1 . 23350E-13 
1 . 55065E-14 
6 . 91224E-16 
2 . 45091E-16 
3 . 08 14  OE-17 
1 . 44896E-18 
4 . 59591E-19 


.2218786602 
.0286304315 
. 0041374068 
.0006441605 
.0001055932 
1 . 790519E-5 
3 . 101577E-6 
5 . 443513E-7 
3 . 053709E-7 
9 . 630086E-8 
1 . 711839E-8 
3 . 051736E-9 
5. A4982E-10 
5 . 48906E— 11 
1 . 74272E-11 
3 . 11853E-12 
4 . 22771E-14 
1 . 00820E-14 
5 . 73400E— 16 
7 . 77566E-18 
1 . 78124E-18 
1. 03795E-19 
1 . 44926E-21 
3 . 41426E-22 


Key: 

Temp  =  Temperature  in  Kelvin 

KA-1  =  Inverse  Temperature  *  1000  in  1/Kelvin 
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Table  18.  Vacancy  Formation  Energy  (Ef) 


Temp (K) 

T (KA-1) 

CuTi2 (eV) 

CuTi(eV) 

FeTi(eV) 

NiTi (eV) 

10000 

.  1 

5000 

.2 

1.03681 

1.17275 

1.31769 

1.76454 

3333.33 

.3 

1.01454 

1.16207 

1.25024 

1.66694 

2500 

.4 

1.00285 

1.15212 

1.20358 

1.60272 

2000 

.5 

.99205 

1.14300 

1.16865 

1.55832 

1666.66 

.6 

.97980 

1.13474 

1.14166 

1.52915 

1428.57 

.7 

.96729 

1.12738 

1.12139 

1.51078 

1250 

.8 

.95551 

1.12093 

1.10662 

1.49948 

1200 

.833333 

.94826 

1.11709 

1.09907 

1.49442 

1111.11 

.9 

.94325 

1.11448 

1.09459 

1.49173 

1000 

1 

.93563 

1.11057 

1.08868 

1.48850 

909.091 

1.1 

.92763 

1.10651 

1.08354 

1.48601 

833.333 

1.2 

.92079 

1.10309 

1.07998 

1.48452 

750 

1.33333 

.91416 

1.09979 

1.07722 

1.48353 

714.286 

1.4 

.90938 

1.09743 

1.07561 

1.48302 

666.667 

1.5 

.90604 

1.09578 

1.07468 

1.48276 

571.429 

1.75 

.90054 

1.09305 

1.07351 

1.48248 

545.454 

1.83333 

.89652 

1.09102 

1.07285 

1.48236 

500 

2 

.89431 

1.08986 

1.07260 

1.48232 

444.444 

2 .25 

.89159 

1.08835 

1.07236 

1.48239* 

428.571 

2 .33333 

.89000 

1.08739 

1.07226 

1.52392* 

400 

2.5 

.88914 

1.08680 

1.07222 

1.46977* 

363.636 

2.75 

.88807 

1.08597 

1.07219 

1.47231* 

352.941 

2.83333 

.  88744 

1.08540 

1.07217 

1.49495* 

333.333 

3 

.88710 

1.08502 

1.07217 

307 . 692 

3 .25 

.88668 

1.08486* 

1.05377* 

300 

3 .33333 

.88643 

1.08245* 

1.18740* 

285.714 

3.5 

.88630 

1.08401* 

*  These  values  are  less  accurate  due  to  their  having  vacancy  concentrations  equal  to  or 
less  than  IE- 18,  thus  their  calculation  errors  are  about  equal  to  the  tolerance  values  of 
the  TK  Solver  model  calculation. 
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APPENDIX  D: 


Stable  Defect  Configuration  Plots 


The  following  figures  are  output  plots  from  the  RUNDYN 
code  for  each  of  the  most  stable  defect  configurations  plus 
example  plots  of  vacancy  and  migration  configurations. 
Figures  12  and  14  are  minimized  vacancy  configurations,  but 
are  not  the  most  stable  defect  configurations.  The  most 
stable  configurations  for  these  lattices  are  Figures  13  and 
15  respectively.  Figures  18  through  22  are  examples  of  Fe 
vacancy  migrations.  The  Ni  vacancy  migration  plots  were 
similar  in  all  respects  except  the  configuration  in  Figure 
21  was  stable  for  NiTi.  The  remaining  Figures  23  through 
38  are  all  of  the  most  stable  interstitial  configurations 
for  either  NiTi  or  FeTi.  In  all  cases  the  Ti  atoms  are 
symbolized  with  circles  with  horizontal  lines  and  the  Ni  or 
Fe  atoms  are  symbolized  with  circles  with  diagonal  lines. 
Except  where  noted  all  of  these  figures  have  unit  direction 
views  of  x=l,  y=0.1,  and  z=0.1.  All  of  the  figures  in  this 
Appendix  were  reproduced  from  the  AFIT  mainframe  system. 

The  rest  of  this  thesis  was  written  using  the  Lotus  soft¬ 
ware  programs  Manuscript  release  2.1  and  Freelance  release 
3.01. 
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Figure  18.  FeTi  Lattice  with  Fe  Vacancy  Direct  Migration 
Atom  #569 
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Fgure  19.  FeTi  Lattice  with  Fe  Vacancy  Migration  Sega 


#1  of  Atom  #440 


Figure  21.  FeTi  Lattice  with  Fe  Vacancy  Migration  Sequence 
#5  of  Atom  #569  is  Unstable  with  Atom  #568  Also  Moving 


Thu  S e f.  2  0  1  2  :  0  ( .  :  3  4  1  9  9  0  :  t  i  b  6  ft  b  +  1  .  0  .  r  e s 


Figure  22.  FeTi  Lattice  with  Fe  Vacancy  Migration  Sequence 
#6  of  Atom  #568 
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Figure  23.  NiTi  Ti  Interstitial  Configuration,  Trajecto 
ries  from  Figure  27 


Thu  Oct.  11 


:  2  2  :  4 1  19  9  0 


Riti  inti 


.  test.4  .  res 


Figure  24.  NiTi  Ti  Interstitial  Configuration  <100>  View 


n  oct  15  1-5:15:00  1  99  0  :  niti_intl.test1.rcs 

Figure  25.  NiTi  Ti  Interstitial  Configuration  <010>  View 


Figure  26. 


NiTi  Ti  Interstitial  Configuration  <001>  View 


Figure  34.  NiTi  Ni  Interstitial  Configuration  <001>  View 


Oct  1C  o:  :  14  : 50  1  990  :  fc.intl . test  3 . res 

Figure  37.  FeTi  Fe  Interstitial  Configuration  <010>  View 


Figure  38. 


FeTi  Fe  Interstitial  Configuration  <001>  View 
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Amorphous 

crystalline  transition,  35 
Antisite 

binding  energy,  18 
code  input  values,  48 
defect  p^opej. ties ,  14 
energy,  35 
formation  energy,  12 
pair  energy,  vii,  35 
structure,  12 
Arhennius  plot,  23 
Atomistic  simulation,  2,  4, 
6,  23 

techniques,  5 

B2-type  Structure 
FeTi,  7 
NiTi,  7 

Binding  energy,  vii,  6,  17 
Boltzmann 
constant,  23 
Boundary 

code  input  conditions,  46 

Chemical  disorder,  14,  35 
Critical  temperature 
Monte  Carlo,  31 

Defect  configurations 
stable  output  plots,  69 
Dumbbell 

Ni-Ni  and  Fe-Fe,  31 
DYNAMO  code,  vii,  5,  42 

Error 

assessment,  10 
TK  Solver  Plus,  61,  68 

FeTi 

iron  titanium,  3 
melting  temperature,  14 
Ti  interstitial,  32 
Fletcher-Powcll 
method,  42 

minimization  method,  5 
minimizer,  48 
Formation  energy 
Fe  antisite  pair,  13 
FeTi , NiTi , CuTi , CuTi2 ,  68 
interstitial,  38 


Ni  antisite,  13 
vacancy,  24,  38,  60 
Frenkel  pair 
creation,  48 

Global  minimum 
configuration,  8 

Input  file 
RUNDYN  code,  44 
Interatomic  potential,  6 
Intermetallic  compounds,  40 
Interstitial 
configurations,  47 
defect  energy,  viii 
defect  properties,  30 
definition,  8 
formation  energy,  38 
migration,  4 
Ni  and  Fe,  34 
stable  configuration,  31 
Ion  microscopy,  23 

Jumping  atom,  43 

Lattice  generation 
code  input,  46 
Lattice  parameter 
FeTi,  7 
NiTi,  7 
Lotus  123 
spreadsheet,  49 
Lotus  software 
Freelance,  69 
Lotus  123,  49 
Manuscript,  69 

Metastable 
configuration,  8 
Migration 
code  input,  47 
Migration  energy,  19 
barrier,  15,  17,  21,  43 
Ni  or  Fe  atoms,  15 
Molecular  dynamics,  5 
Molecular  statics,  5 
Monte  Carlo,  8 
code  input,  48 
critical  temperature,  30 
technique,  5,  8,  30 
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NiTi 

melting  temperature,  14 
nickel  titanium,  3 
Ti  interstitial,  3*.’ 

Output  file,  45 
Output  plots,  69 

Perfect  lattice 
structure,  11 
system  energy,  11 
Point  defect 
definition,  1 
Potential  energy,  5 
code  input,  45 

Quenching 
metal,  23 

Reaction  coordinate,  43 
Restart  file,  45 
RUNDYN  code,  5,  9 
description,  42 

Significant  digits,  10,  49 
Simulated  annealing,  5,  8, 

30 

Softlink,  44 
Source  code,  44 
Split-interstitial,  viii,  37 
Cu-Cu,  39 
Ni-Ni  dumbell,  31 
Ni-Ni  or  Fe-Fe,  41 
Structure 
FeTi,  7 
NiTi,  7 

Ti  Interstitial 
FeTi  configuration,  32 
NiTi  configuration,  32 
Ti  Vacancy 
configuration,  16 
energy,  15 
TK  Solver  Plus 
iterative  solver,  61 
software,  27 
software  code,  56 
Tolerance,  10 
RUNDYN  code,  48,  49 
TK  Solver  Plus,  61,  62,  68 
Transition 

crystal 1 ine-to-amorphous , 

35 


Vacancy 

concentration,  23 
concentrations,  67 
configurations,  35 
CuTi  concentration,  65 
CuTi2  concentration,  66 
defect  energy,  vii 
defect  properties,  19 
definition,  8 
FeTi  concentration,  64 
formation  energy,  6,  36, 
60,  86 

migration,  vii,  36 
migration  energy,  22 
migration  sequence,  19 
Ni  or  Fe,  18 
NiTi  concentration,  63 
Ti  vacancy,  15 
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